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ABSTRACT 



A computer model which simulates Naval Gunfire Support 
of Small Independent Action Forces (SIAF) in an insurgent- 
contested area is presented. The mechanics of the model are 
described in sufficient detail ^o allow understanding by any 
user with a basic understanding of Naval Gunfire Support and 
computer simulation techniques. A detailed guide to pre- 
paring the input data deck required for simulation of a user- 
defined problem is included. 

The problem definition and results of a partial fea- 
sibility study of the Seabased Tactical Deterrent Force 
(STDF) concept using the simulation model are presented as 
a means of providing some insight into the feasibility of 
the STDF concept and as an example of the use of the simula- 



tion model . 
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I. INTRODUCTION 



In these days of reduced U.S. Military manpower and 
hardware, a great deal of attention is being focused on 
developing the optimal military posture for future opera- 
tions. One of the main factors being considered in devel- 
oping this posture is the United States' announced desire 
to reduce the number of U.S. troops and troop support fa- 
cilities on foreign soil while still maintaining the ability 
to respond with an adequate amount of military power to meet 
commitments to allied countries. Thus, a foreseeable pos- 
sibility is that the U.S. might be called upon for emergency 
military aid to a country in which it has no land-based 
troops, aircraft, or artillery, as was the case in the 
Lebanon and Dominican Republic crises of the recent past. 

One of the concepts of future military posture which is 
under consideration is that of the Seabased Tactical Deter- 
rent Force (STDF) . Integral to the STDF concept is the con- 
cept of the Small Independent Action Force (SIAF) . Each of 
these concepts is discussed below. 

A. SEABASED TACTICAL DETERRENT FORCE 

The STDF would be required to have a multitude of cap- 
abilities. Without defining all of the capabilities re- 
quired, it can be summarily stated that the STDF would have 
to combine the existing capabilities of an Attack Carrier 
Striking Force, an Amphibious Task Force, a Naval Gunfire 
Support Force and an Underway Replenishment and Logistics 
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Force, while providing for its own defense as well. One of 
the most important requirements of the STDF, and the one 
which emphatically distinguishes it from other alternative 
concepts, would be the ability to strategically interpose 
combat forces in most areas in which conflicts might arise 
which would seriously jeopardize vital U.S. interests. This 
is not a new requirement, but it would enjoy renewed and 
expanded importance under the STDF concept. The convention- 
al approach to exercising this capability which has been 
used many times in the past has been to land large numbers 
of troops ashore by an amphibious landing to conduct massive 
operations to achieve a military objective. But with 
changing U.S. foreign policy, such massive operations ap- 
pear less likely. Also, our experience in Southeast Asia 
has inspired new technology and concepts of amphibious war- 
fare. In the face of an insurgent force employing guerrilla 
warfare or similar tactics, massive amphibious operations 
are rarely practical or productive. The effectiveness of 
small reconnaissance and combat patrols, such as the Small 
Independent Action Force (SIAF) patrols, in seeking out the 
insurgent forces has often been demonstrated. Thus, the 
insertion and support of SIAF patrols becomes a capability 
which the STDF must have. 



B. SMALL INDEPENDENT ACTION FORCE (SIAF) PATROLS 

At this point in the discussion it is presumed that it 
has been determined that the STDF concept is a viable means 
of interposing U.S. military power, and the ability to insert 
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and support SIAF patrols is organic to the STDF. The fol- 
lowing question then arises; "What type of operations by 
the SIAF patrols would be most effective in a limited \'iar 
situation if one arises?" Some high-ranking military of- 
ficers have expressed the opinion that it might be tactic- 
ally advantageous in many areas to saturate an area of 
limited size with SIAF patrols. These patrols would be 
dispersed throughout the area with the objective to either: 

1. Reconnoiter an area where enemy strength and degree 
of control is virtually unknown, or 

2. Seek out known or suspected enemy forces. 

In either case, each individual patrol should typically 
be of short duration (a few days) , small in numbers (4 to 18 
men), dependent only on the STDF for logistics and external 
fire support, self-sustaining while actually on patrol, and 
capable of rapid and reliable communication with the STDF. 

The primary mission envisioned for the SIAF patrols would 
be one which is essentially passive. The patrols would 
seek out enemy personnel, base camps, support facilities, 
etc. , but rather than actually engaging in combat with such 
targets, the patrols would call in air strikes or naval gun- 
fire from the STDF. Thus, the mission of the patrols is 
primarily reconnaissance/forward observer. 

C. OBJECTIVES OF THIS THESIS 

In this paper, a small segment of the overall operations 
analysis effort required to investigate the STDF and SIAF 
concepts is isolated and conducted. Specifically, the subject 
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under study is Naval Gunfire Support (NGFS) , originating 
from the STDF, in support of SIAF's ashore in an insurgent- 
contested area. The objective of the paper is twofold: 



1 . Development of Computer Simulation Model 

A computer simulation model has been developed 
which simulates the actions of the STDF in providing NGFS 
to SIAF's ashore. The model is designed primarily for in- 
vestigating the queueing aspects of providing NGFS. The 
main focus of attention is the time interval between the 
time a SIAF calls for a mission to be fired and the time 
a ship commences firing the mission. No ascertainment is 
made of the actual results achieved I'/hen the mission is 
fired. The model is very flexible in design, allowing the 
suTiUlatron of a wxde varxcty of sxtuatxons dcfxncd by user- 
specified inputs. 

2 . P artial STDF Feasibility Study 

After the simulation model had been developed, a 
typical situation in an insurgent-contested area was hy- 
pothesized, and the model ivas exercised to determine the 
optimum number and mix of various types of NGFS ships which 
would be required to adequately support the SIAF's ashore. 
Since situations in different areas certainly vary greatly 
and may be more or less severe than the hypothesized situ- 
ation, and since the simulation is restricted solely to NGFS 
and no other organic capabilities of the STDF (such as air 
strikes) , the results of the exercise cannot be considered 
to provide any definite conclusion concerning the overall 
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feasibility of the STDF concept. However, the exercise 
provides an example of the use of the simulation model for 
other users to follow when a specific real-life situation 
is under investigation, and it also provides some insight 
into the feasibility of a STDF providing NGFS in support 
of SIAF’s ashore in an insurgent- contested area. 

D. ORGANIZATION OF THIS THESIS 

The development and use of the model are discussed in 
this thesis. In the next section, the systems to be simu- 
lated are described to the degree of detail necessary to 
construct the simulation model. Then, the mechanics of the 
model are discussed in sufficient detail to allow any user 
with a basic knowledge of NGFS and simulation techniques 
to understand how the computei' program works. Next, the 
method and results of a partial STDF feasibility study are 
presented and conclusions based on the results are dis- 
cussed. Finally, refinements and extensions of the model 
are discussed. 
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II. THE SYSTEM 



The primary system which is simulated can be identified 
as Naval Gunfire Support . In order to construct a realistic 
and useful simulation model, at least two other systems must 
be included in the simulation model. These are refueling 
and rearming of ships at sea and the conduct of SIAF patrols 
ashore. Each of these systems is described in this section 
of the paper. 

A. NAVAL GUNFIRE SUPPORT 
1 . NGFS Ship Types 

Any U.S. Naval ship v^?ith installed guns is considered 
capable of providing some degree of gunfire support to troops 
ashore. However, from among the many types of U.S. Naval 
ships ^^;ith installed guns, the only ones which can be con- 
sidered to have NGFS as one of their primary duties are 
cruisers (CA, CL, CLG) , destroyers (DD, DL, DDG , DLG) , and 
battleships (BB) . In recent years, certain types of ships, 
such as the Inshore Fire Support Ship (IFS) and the Medium 
Rocket Launching Ship (LSMR) , have been built or converted 
with the specific primary mission of providing NGFS. Hov^^- 
ever, these ships normally do not have adequate cruising 
ranges and speeds to be considered as permanent elements of 
a STDF. There are no battleships in commission at the 
present time. Therefore, the only types of NGFS ships con- 
sidered in the model are cruisers and destroyers. 
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2 . 



Installed Guns 



All cruisers and destroyers have more than one 
type of gun installed, but each has one type which is con- 
sidered to be it's "main battery" for NGFS. For example, 
a heavy cruiser (CA) normally has 6 8"/55 guns, 10 5"/38 
guns, and 8 3"/50 guns. However, the 8"/55 guns are con- 
sidered to be the main battery and would normally be the 
only ones used for NGFS, while the other guns are primarily 
used as a backup and for self-defense in surface-to-air 
and surface-to-surface engagements. In the model, each 
ship is considered to use only its main battery for NGFS. 
The types of ships considered in the model and the main 
battery gun type of each [Reference 1] are listed below: 



Cruis 


ers 


Destroyers 




CA 


8"/55 


DD (Fram I/II) 


5"/38 


CL 


6"/47 


DD (931 Class) 


5"/54 


CLG 


6"/47 


DDG 


5"/54 






DLG 


5"/54 






DL 


5"/54 



Throughout this paper, a CA is referred to as an 
8"/55 ship, a CL or CLG as a 6"/47 ship, etc. 

3 . Firing Procedures 

Procedures for calling for and providing NGFS vary 
with the specific situation and local instructions. In some 
situations ships might simply proceed to a firing position 
and fire upon a predetermined area without the assistance of 
spotters ashore. This is called harrassment and interdiction 
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(H§I) fire, but this type of fire is not envisioned in the 
STDF/SIAF concept. The type of firing which is most likely 
to occur in the STDF/SIAF situations can be considered in 
almost all cases to adhere to the following chronological 
sequence of events : 

a. Call for Fire 

A SIAF patrol detects a target and radios a 
request for NGFS to the STDF. This is referred to as a 
"call for fire”, and will be abbreviated as CFF throughout 
this paper. 

b. Ship Assigned 

The STDF commander or his delegated subordinate 
determines which ship from among those available, if any are 



available , 


is mo s 


t 3.pp rep ri 3.1^6 “to fire the ini ss ion ap_d. 9.s~ 


signs that 


ship . 


This assignment is based on the following 


factors : 








CD 


Range to the target. 




(2) 


Type of the target. 




(3) 


Distance of each ship from the firing posi 



tion . 

c. Assigned Ship Proceeds to Firing Position 

The firing position fbr each mission is considered 
to be the point along the "gunline” which is closest to the 
target (the perpendicular distance from the gunline to the 
target) . The gunline is a fictitious line in the ocean run- 
ning roughly parallel to the shoreline at a specified dis- 
tance from the shore. 
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d. Assigned Ship Sets Up for Mission 

After the assigned ship has reached the firing 
position, it fixes its navigational position and readies 
its gun and then reports to the SIAF that it is ready to 
fire . 

e. Fire the Mission 

The assigned ship commences firing on a signal 
from the SIAF. After the first few rounds have been fired, 
the SIAF observes the fall of shot and radios instructions 
to the ship to bring fire onto the target. This is called 
"spotting”. l\nien fire has been adjusted onto the target, 
the ship continues to fire at a rapid rate until the SIAF 
observes that the target has been destroyed or otherwise 
effectively encountered. The ship then ceases fire and re- 
turns to a "ready” status. 

4 . Ship Operations When Not Firing 

Ships assigned to NGFS duties but not actually en- 
gaged in firing would normally remain in the vicinity of the 
gunline in order to be able to respond rapidly when assigned 
to fire a mission. The specific area in which non-firing 
ships would operate would be specified by the STDF commander. 
A frequently used procedure is for non-firing ships to oper- 
ate on or inside of a race-track type waiting pattern at a 
distance from the gunline determined by considering the 
tradeoffs between rapid response time and the need to keep 
non- firing ships out of range of suspected enemy shore bat- 
teries . 
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B. REFUELING AND REARMING 

A definite limitation on the availability of NGFS ships 
is imposed by fuel and ammunition capacities of each ship. 
Under normal steaming conditions, a destroyer with a full 
load of fuel can steam for 3 to 4 days before it will be re- 
quired to refuel, while a cruiser can steam 5 to 6 days. 

This assumes that refueling is required when the amount of 
fuel on board reaches 50% of ships capacity, a fairly com- 
mon criterion. To determine when rearming is required, the 
STDF commander would normally specify a minimum on-board 
ammunition level to which ships would be allowed to fall be- 
fore rearming. This level would be a function of the type 
of ship, expected rate of fire, and number of other ships 
present. 

It is possible that refueling and rearming requirements 
would not impose a limitation on NGFS provided. This v\fould 
be the case when other NGFS ships are available to tempor- 
arily replace ships which are required to refuel or rearm, 
or when the duration of the operation is short enough that 
rearming and refueling should not be necessary. 

C. SIAF PATROL PROCEDURES 

A very detailed computer simulation model of SIAF patrol 
actions has been developed by TRW, Inc. of Los Angeles and 
is described in detail in Reference 2. This model can be 
used for developing optimal patrol procedures and for deter- 
mining expected patrol results for a wide variety of friendly 
and enemy situations. In the model described in this thesis. 
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no attempt has been made to simulate specific patrol actions 
with the exception of patrol insertion, calls made to the 
STDF for NGFS, and patrol extraction. Under the STDF con- 
cept, a SIAF patrol would typically adhere to the following 
chronological sequence of major events: 

1 . Insertion 

The patrol is inserted into its assigned subarea of 
the overall objective area. This insertion would normally 
be accomplished by helicopter from an LHA assigned to the 
STDF, but could conceivably be accomplished by other means 
such as paradrop or amphibious landing. Order and time of 
insertion into the various subareas would be chosen such as 
to yield the highest probability of achieving a tactical ad- 
vantage by surprising and confusing the enemy. 

2 . Patrol Maneuvers 

The SIAF patrol maneuvers within its assigned sub- 
area attempting to locate the enemy and his support facil- 
ities while remaining undetected by the enemy. When targets 
are located the patrol calls for NGFS or air strikes, if 
available, rather than actually engaging in combat with the 
enemy. The patrol would continue to maneuver within its as- 
signed subarea until extracted forgone of the three reasons 
indicated below. 

3 . Extraction 

The patrol is extracted when one of the following 

occurs : 

a. A predetermined extraction time is reached. 
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b. The patrol leader considers that optimal area 
coverage has been achieved. 

c. Extraction is required because the patrol's 
effectiveness has been significantly reduced by personnel 
casualties . 

Extraction would normally be accomplished by heli- 
copter. 

D. THE ENEMY SITUATION 

Obviously, the enemy situation is a factor which plays 
a large role in any study of NGFS and SIAF patrols. No 
study could ever begin without a reasonably accurate estima- 
tion of the size of the enemy force, its locations, its cap- 
abilities, and its expected actions. These enemy factors 
and others are considered in great detail in the TRW SIAF 
model. They are considered in the model described in this 
paper in the following ways: 

1 . Subarea Division 

Division of the objective area into patrol sub- 
areas is based on terrain features, expected enemy locations, 
and expected number and capabilities of enemy in each sub- 
area. 

2 . Call - for- Fire Rate 

The real cornerstone on which the model depends for 
accurate estimates of the queueing aspects of NGFS is the 
expected rate of calls for fire in each subarea. This rate 
is, of course, primarily dependent on the enemy density in 
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each subarea. For a specific situation, the CFF rate for 
a particular subarea can be closely estimated by the TRW 
SIAF model or can be predicted based on available intelli- 
gence . 

E. ASSUMPTIONS 

There are several general assumptions which are neces- 
sary in modeling NGFS within the scope of this paper. These 
are discussed below. V/hile each of these detracts from the 
realism of the model, it is considered that none of them 
significantly detracts from the accuracy of queueing re- 
sults. Other more specific assumptions will be discussed 
as they occur in the description of the model in Section 
III . 

1. SIAF patrols remain within their assigned subarea 
throughout the time interval between insertion and extrac- 
tion, and only detect targets in their assigned subarea. 

2. The enemy situation at the start of the operation 
in each subarea remains the same throughout the operation. 

3. The enemy does not launch attacks against NGFS 
ships, or if he does the attacks are not successful enough 
to significantly reduce the operational capability of the 
ship . 

4. SIAF patrols do not come under enemy attack, or if 
they do the attack is such that emergency NGFS is not re- 
quired and the patrol's capability is not seriously reduced. 

5. A SIAF patrol remains stationary at the same point 
from the time it has called for a mission until the mission 
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is complete. This leads to the further assumption that 

t 

patrols spot only one target at a time, i.e., while one 
mission is in progress the patrol does not detect any 
other targets. 

6. IVhen a SIAF patrol is extracted from its assigned 
subarea, another patrol is not inserted into the area to 
replace the extracted patrol. 

7. SIAF patrols do not sleep; at least, no element of 
the model specifically simulates a period of time during 
which patrols are sleeping. However, it can be considered 
that unusually long interarrival times between missions in 
a particular subarea could be caused in some instances by 

a period of time during which the patrol in that subarea is 
sleeping . 
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III. THE MODEL 



This section describes the model which has been devel- 
oped to simulate the systems described in Section II. The 
simulation program is written in FORTRAN IV. It is not the 
intent of this section to describe the entire program step- 
by-step; but rather to describe the mechanics of the program 
in sufficient detail to allow any user with a basic knowledge 
of NGFS and computer simulation techniques to understand what 
is being simulated and to be able to use the model according 
to his own input parameters, or to adapt certain sections of 
the program, if desired. A generalized flow chart of the 
program is contained in Appendix A and the complete program 
listing is contained in Appendix B. Numerous comment cards 
have been included in the program listing to facilitate un- 
derstanding . 



A. PRELIMINARY DISCUSSION 

1 . Variable Names 

Mnemonic variable names are used throughout the pro- 
gram to facilitate association between the names and the 
quantities which they represent. 

2 . Measures 

All times in the program are measured in whole 
minutes and distances are measured to the nearest whole yard, 
unless otherwise indicated. In view of these measures, all 
variables beginning with A through W are declared to be 4- 
byte integers. Whenever a calculation of an integer-valued 
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time quantity is performed, 0.5 is added to the result, re- 
sulting in a truncated integer accurate within ±.5 minutes. 

3 . Input Data-General 

Almost all numerical quantities used in program cal- 
culations are read in at the beginning of the program from a 
data deck prepared by the user. The definition and use of 
each input variable is discussed as it appears in the model 
description which follo^^/■s. A detailed guide to assembling 
the data deck is included in Section II I. J. 

4 . Random Numbers 

There are 14 points in the program at which random 
numbers are used for Monte Carlo-type comparisons or calcula- 
tions. Random number streams at each point are independent 
of all other points in the program. The random number gen- 
erator is initialized by specifying any 14 odd five-digit 
numbers as input variables IRl through IR14. KR is set equal 
to 65539 at the beginning of the program. Random numbers 
are then produced by the following calculation at each point 
a random .number is required (IRl is used for example): 

IRl = IRl * KR 

XRN = 0.5 + FLOAT(IRl) * 2.328306E-10 

XRN will be a floating point random number uniformly 
distributed between 0.0 and 1.0 [Reference 3]. 

5 . Ship Identification 

Throughout the program ships are referred to by 
their type (T) and ship number (S) . The four types of ship 
and corresponding values of T are: 
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1 = 8"/55 ship (cruiser) 

2 = 6”/47 ship (cruiser) 

3 = 5''/54 ship (destroyer) 

4 = 5"/38 ship (destroyer). 

A maximum number of 9 ships of each type are pro- 
vided for, or a maximum of 36 ships. The number of each 
type of ship ’’present in the problem” is specified by the 
input variables NR855 , NR647 , NR554, and NR538. ’’Present 
in the problem” refers to the number of ships present at 
the start of the problem. The names of quantities which 
are related to a particular ship are subscripted by the 
ship number and type in that order, i.e., UTIL (5,2) is the 
cumulative total minutes that 6”/47 ship tfS has been uti- 
lized for NGFS. 

6 . Uniform Distributions 

There are several points in the program where ob- 
servations of a random variable from a uniform distribution 
are made. In each case, the observations are generated by 
use of the following formula: 

RV = MIN + XRN * (MAX - MIN) 

where RV is the random variable being observed, XRN is a ran- 
dom number between 0.0 and 1.0, and MIN and MAX are the pa- 
rameters of the uniform distribution. 

7 . Number of Iterations 

The user specifies the number of iterations of the 
problem which he desires to use to obtain averaged results. 
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This number is specified as the input variable NRIT. In 
several trial runs of the program, NRIT = 10 was sho\m to 
be sufficient for obtaining reliable results with a rela- 
tively narrow 95% confidence interval. However, increasing 
NRIT will reduce the width of the confidence interval even 
further. 

8 . Specification of Card Reader and Printer Logical 

Unit Numbers 

All READ and WRITE statements in the program use th 
values of the variables R and W as the card reader and prin- 
ter logical unit numbers, respectively. As written, the 
program specifies R = 5 and W = 6 at the beginning by the 
statement DATA R,W/5,6/. If the card reader and printer 
logical unit numbers at a computer facility on which the 
program is to be run are not 5 and 6, respectively, then the 
above DATA statement must be replaced by one which specifies 
the appropriate logical unit numbers. 

B. PROBLEM AREA 

The first step in developing data for the model is to 
define the problem area in terms of a coordinate system 
which can be specified to the computer. 

1 . Objective Area CQA) 

The OA is the overall area ashore in which SIAF 
patrols are conducted. The OA is divided into NRAREA patrol 
subareas, where NRAREA is an input variable specifying the 
number of subareas. NRAREA must be bet^^^een 2 and 50. These 
subareas must be rectangular and the boundaries of each sub- 
area must be parallel or perpendicular to the boundaries of 



all other subareas, and the northern and southern boundaries 
of each subarea must be parallel to the gunline. For pur- 
poses of uniformity it is considered that, all subarea bound- 
aries run north to south and east to west. The western edge 
of the westernmost subarea is considered to be the y-axis of 
an x-y coordinate system (see Section III.B.3 for exception) 
The distance of the eastern and western edges of each sub- 
area from the y-axis are the input variables EAST(NR) and 
WEST(NR), where NR is the subarea number. Subareas may be 
numbered in any way desired by the user. 

2 . Gunline 

The gunline is a fictitious straight line in the 
ocean, parallel to a straight line oriented east-west which 
approximates the shore line. See Section II.A.3.C for a 
definition of the gunline. The gunline is considered to be 
the x-axis of the x-y coordinate system. The distance of 
the northern and southern boundaries of each subarea from 
the gunline are the y-coordinates of the subarea, specified 
by the input variables NORTH (NR) and SOUTH (NR). The y- co - 
ordinate of the northern edge of the northernmost subarea 
must be greater than the maximum range of the longest gun 
present . 

3 . Waiting Pattern 

The area in which ships operate while not assigned 
to a mission is defined by the user by appropriate specifica 
tion of values of the following input variables: 

WTDIS - Perpendicular distance of the northern edge 
of the waiting pattern from the gunline 
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NSWAIT - The north-south length across the waiting 

pattern 

WTNW and WTNE - The perpendicular distance of the 
western and eastern edges, respectively, of the waiting pat- 
tern from the y-axis. (NOTE: If the western edge of the 

waiting pattern is west of the western edge of the OA, then 
the western edge of the waiting pattern should be used as 
the y-axis (WTNW = 0)). 

4 . Justification of Problem Area Definition 

It is realized that the scheme of problem area 
definition described above does not exactly conform to most 
real-life situations. For example, area boundaries would 
often conform to rivers or treelines and not necessarily be 
straight lines. However, straight line approximations such 
as used in this model should not significantly effect the 
queueing results if input variables are properly specified. 

Figure 1 illustrates a typical problem area defined 
in the coordinate system described in this section. 



C. EVENT CHAIN AND .MISSION SCHEDULE 

The two principle mechanisms by which data is kept track 
of and program control is exercised are the "Event Chain" and 
the "Mission Schedule". 

1 . Event Chain 

The simulation uses the "next-event" method of time 
advance. At the beginning of the problem the statements 
K = 1 and CLOCK = EVENT(K) sets the problem timer to the time 
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Example of Area Coordinates: 



4 



NORTH(l) = 20000 
SOUTH(l) = 13000 
WEST(l) = 0 
EAST(l) = 8000 



True Problem 
North North 



Figure 1. Typical Problem Area Map. 



32 




4 



of the first event. This event is then directed through its 
appropriate path of logic, l^dien it has been completely pro- 
cessed the value of K is increased by 1 and the time is ad- 
vanced to the time of the next event. This process is 
continued until an event \>:ith a time of 9999999 is the next 
event and the problem is terminated. Such an event is 
created when the generated time of the next call for fire 
in a subarea is greater than the time of extraction of the 
patrol from that subarea. 

a. Event Types 

There are four types of events considered in 
the program. These are indicated by the value of the vari- 
able EVTYP(K) as follows: 

1; A call for fire from a SIAF patrol 
2: A NGFS mission completion 

3: A scheduled refueling of a NGFS ship 

4: A completion of refueling, rearming, or 

combined refueling/ rearming of an NGFS ship. 

b. Event Parameters 

The Event Chain actually consists of four param- 
eters associated with each event. These are stored in vector 
arrays with variable names as follows: 

EVENT(K) : Time of event (minutes from zero 

problem time) 

EVTYP(K): Type of K^h event (values lAW Section 

Ill.C.l.a) 

EVTMSN(K): If EVTYP(K) = 1 or 2 , then EVTMSN(K) 
is the mission number of the K^h event. Mission numbers are 



assigned in chronological order (see next subsection). If 
EVTYP(K) = 3 or 4, then EVTMSN(K) = 0. 

REPSHP(K): If EVTYP(K) = 1 or 2 , then REPSHP(K) 

= 0. If EVTYP(K) = 3 or 4, then REPSHP(K) is a two digit 
number representing the number and type of the refueling 
and/or rearming ship, i.e., REPSHP(K) = 23 m.eans 5"/54 ship 
#2 is either scheduled to refuel or has completed refueling 
and/or rearming. 

c. Subroutine VENTCH 

Each time a new event is generated in the pro- 
gram subroutine VENTCH is called. The arguments of VENTCH 
are the four parameters discussed above. S/R VENTCH inserts 
the new event into its proper chronological order on the 
event chain, and increases the event number of all suceeding 
events by 1 . 

2 . Mission Schedule 

The mission schedule represents a chronological 
listing of all NGFS missions which have been generated dur- 
ing the problem. 

a. Mission Generation 

A call for fire represents a mission arrival. 

In view of assumption 5 in Section* I I . E . , mission interarri- 
val times are represented by the time interval between the 
time the last mission in a particular subarea was completed 
and the time of the next mission arrival in that subarea. 
Mission interarrival times in each subarea are assumed to be 
exponentially distributed with a rate which is indirectly 



specified by the user in the following manner: For each sub- 

area, the minimum and maximum expected patrol duration (num- 
ber of days from insertion time to the nearest tenth) are 
specified by the user as the input variables XSHORT(NR) and 
XLONG(NR). The program calculates expected patrol duration 
as 

XDURAT(NR) = (XSHORT(NR) + XLONG (NR) ) /2 . 

For each subarea, the user also specifies the expected number 
of calls for firing during a patrol (i.e., during the time 
XDURAT(NR)), as the input variable XCFF(NR) . Then the pro- 
gram calculates the call- for-f ire rate for each subarea as 

XCALRT(NR) = XCFF (NR) / (XDURAT (NR) * 1440) 

where the factor 1440 converts days to minutes. The arrival 
rate then has units of calls/minute. 

Exponentially distributed interarrival times 
are generated as follows [Reference 4] ; 

The cumulative distribution function for the 
exponential distribution is 

F(x) = 1 - e‘^^ 

where x is an interarrival time and X is the arrival rate. 
Since F (x) is a continuous probability between 0.0 and 1.0 
it can be replaced by a random number between 0.0 and 1.0, 
say r. Then 

r = 1 - e‘^^ =5>x = - (l/X)ln(l-r) 
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or equivalently 



X = - (1/ A) ln(r) . 

Thus, -ALOG (XRN) /XCALRT (NR) is an exponentially distributed 
interarrival time of a call for fire in subarea NR. There 
are two points in the program at which calls for fire are 
generated . 

(1) First Mission . At the start of the problem, 
the time of the first mission to be called for in each sub- 
area is calculated as 

FIRST (NR) = INSERT (NR) - ALOG (XRN) /XCALRT (NR) +0.5 

where INSERT(NR) is the time of patrol insertion into subarea 
NR in minutes from zero problem time. INSERT(NR) is an in- 
put variable specified by the user. Zero problem time is 
defined as the number of minutes after midnight on the first 
day of the problem that the first partrol is landed in its 
assigned subarea. It is specified by the input variable 
HHOUR. The value of HHOUR is used in the program only to 
convert times of calls for fire and mission completion to 
date-time -groups in the Mission History printout. 

These first missions are then placed on 
the mission schedule by setting MSNSKD(I) = FIRST(I) and 
arranging the elements of MSNSKD in chronological order by 
a ''ripplesorf technique. The first missions are then placed 
on the Event Chain with parameters EVENT(K) = MSNSKD(I) , 
EVTYP(K) = 1, EVTMSN(K) = I, and REPSHP(K) =0. 
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(2) Next Missions. Each time a mission is com- 



pleted, the time o£ the next call for fire in that subarea 
is calculated as 

NEXT = CLOCK - ALOG (XRN) /XCALRT(AREA(E;) ) +0.5 

where CLOCK is the present value of the problem time and 
AREA(E) is the subarea in which the mission was just com- 
pleted. This mission is then placed on the mission schedule 
by a routine which is similar to VENTCH but which is not a 
separate subroutine since it is only used at this one point 
in the program. Next missions are also placed on the Event 
Chain with appropriate parameters, 
b. Mission Parameters 

As with the Event Chain, the Mission Schedule 
actually consists of 19 parameters which are assigned to 
each mission. Wien a mission is first generated it is only 
assigned 2 parameters, MSNSKD(I) and AREA(I) . The other 17 
parameters are assigned at various points in the program 
after CLOCK time has reached the time of the call for fire 
for the mission. A discussion of each of the vector ar- 

rays in which these 19 parameters are contained and the way 
in which each is assigned follows;^ 

(1) MSNSKD(I) : Time of call for fire for the 

mission (minutes from zero problem time) 

(2) AREA(I) : Number of the subarea in which 

the mission occurs 

(3) RANGE (I) : Range from the firing position 

to the target of the mission (the target's y- coordinate) 
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(4) EWTAR(I): East-west coordinate of the tar- 
get of the mission (the target's x- coordinate) 

Although target densities vary in different 
subareas, targets are considered to be uniformly distributed 
within each subarea. Therefore the target coordinates are 
calculated as uniformly distributed random variables between 
SOUTH(A) and NORTH(A) for RANGE(I), and WEST(A) and EAST(A) 
for EWTAR(I), where A = AREA(I). 

(5) TYPE(I) : Type of target of the I^h mission 

There are four target types considered in the 

problem. These are; 

1: Many personnel (20 or more) or personnel 

widely dispersed. 

2: Few pei'sonnel not widely dispersed. 

3: Heavy material objects (large bunkers, tanks, 

heavy artillery, etc.). 

4: Light material objects (small bunkers, jeeps, 

etc . ) . 

The probability of each type of target (i.e., 
the estimated percentage of all targets in the subarea which 
are of each type) for each subarea is specified by the user 
by the input variables XPROBl(NR), XPROB2 (NR) , XPR0B3(NR), 
and XPROB4(NR). The sum of these four variables must equal 
1.0 for each subarea. When a call for fire occurs, a random 
number is generated and compared with the cumulative target 
type probability distribution for the area in which the mis- 
sion occurs to determine the target type. 



(6) MISTYP(I) : Mission type of the ith mission. 

This type refers to the types of ships which may be assigned 
to fire the mission, as will be discussed in Section III.D. 

(7) TYPSHP(I) : Type of ship assigned to fire 

the mission. 

(8) SHIPNRCn •• Number of the ship assigned to 
fire the Ith mission. 

The ways in which parameters (7) and (8) are 
assigned are indicated in Sections III.E. and III.F. 

(9) TRAVEL (I) : The time interval between the 

time a ship is assigned to fire the I^h mission and the time ' 
the ship arrives at the firing position. 

(10) SETUP (I) : The time interval between the 

time the assigned ship arrives at the firing position and the 
time the first round is fired. 

(11) FI RE (I) : The time interval between the 

time the first round is fired and the time the last round 
is fired. 

(12) SERVC(I) : Total service time of the 

mission. SERVC(I) = TRAVEL(I) + SETUP(I) + FIRE(I) . 

(13) AMMEXP (I) : Number of rounds of ammunition 

expended in firing the I^h mission. 

Parameters (9) through (13) are all computed by 
subroutine SERVIC, which will be discussed in Section III.G. 

(14) FSTSHT(I) : The time between the call for 

fire and the time the first round was fired for a mission 
which required no queue time. FSTSHT(I) = TRAVEL(I) + SETUP(I). 
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(15) FSTSHO(I): Same as FSTSHT(I) for amission 

which required queue time. FSTSHO(I) = QTIME(I) + TRAVEL(I) 

+ SETUP(I) . 

(16) JOINQ(I): The time the Ith mission joined 

the queue. JOINQ(I) = 0 if the mission required no queue 

time . 

(17) LEAVEQ(I): Time the I^h mission left the 

queue. LEAVEQ (I) = 0 if the I^h mission required no queue time. 

(18) QTIME(I): Amount of time the Ith mission 

was required to wait for assignment of a ship to fire the 
mission (queue time). QTIME(I) = LEAVEQ(I) - JOINQ(I). 

(19) MSNCMP(I): The time at which the 

mission was completed. MSNCMP(I) = MSNSKD(I) + QTIME(I) 

+ SERVC(I). 

D. SHIP ASSIGNMENT DOCTRINE 

1 . Types of Ship Assignment Doctrines 

There are two types of ship assignment doctrine 
available to the user. These refer to the decision process 
which the STDF commander uses to determine which types of 
ships from among those "present in the problem" are appro- 
priate to fire a particular missio^. The type of ship as- 
signment doctrine is specified by the input variable ASNDOC, 
which has values defined as follows: 

1: " Any ship with guns of adequate range may fire ." 
The target range is the only consideration in determining 
which types of ships are appropriate to assign to fire the 
mission. 
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" Best effect on target .” Both target range 



and target type are considered in determining which types 
of ships are appropriate to assign to fire the mission. 

When ASNDOC = 2 is used, the following rules govern 
the determination of appropriate ships: 

1. Heavy material targets (type 3) may only be 
fired at by 8"/55 or 6"/47 ships, unless there are no 8"/55 
or 6”/47 ships present in the problem or there are no 8"/55 
ships present and the target is out of range of a 6"/47 gun. 
The reason for this rule is that 8"/55 and 6"/47 shells have 
a considerably higher explosive po\\rer than 5"/54 or 5"/38 
shells . 

2. Personnel targets (types 1 or 2) may only be 
fired at by 5"/54 or 5"/38 ships and a 5"/54 ship is prefer- 
rable due to its rapid rate of fire, unless there are no 
5"/54 or 5"/38 ships present, or the target is out of range 
of a 5"/54 gun, or there are no 5"/54 ships present and the 
target is out of range of a 5"/38 gun. 

2 . Mission Types 

Regardless of which assignment doctrine is used, 
each mission is assigned a mission type number (MISTYP(I)). 
This number indicates what types of ships, from among those 
present in the problem, are appropriate (eligible) to fire 
the mission. There are 8 possible values of the variable 
MISTYP(I), as follows: 
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MISTYP(I) 



Type Ships Eligible to Fire 



1 



8'755 



2 



8"/55, 6"/47 



3 



5"/54 



6 



4 



5 



7 



5"/54, 5'738 
Any 

8’755, 6’747, 5’754 
8’755, 5’754 



8 



None 



3 . Maximum Gun Ranges 

The maximum range of each type gun in the problem 
is specified by the user as the input variables MAX855, 
MAX647, MAX554, and MAX538. The maximum range of the longest 
type of gun present is also specified by the input variable 
MAXNGF. Due to the construction of the model there are two 
requirements on these quantities: 

'1. MAX855 < MAX554 < MAX647 < MAX538, unless no 
ships of a particular type are present (i.e., if NR554 = 0, 
then MAX554 = 0) . 

2. MAXNGF < MAX {NORTH (NR)} 

NR 

In view of the first requirement, the only type of 
missions which will occur if ASNDOC = 1 are 1,5,6, and 7. 

Any of the types 1 through 7 may occur if ASNDOC =2. A 
type 8 mission will never occur if the second requirement 
is satisfied. If a type 8 mission does occur, the message 

"A TYPE 8 MISSION HAS OCCURRED IN AREA " will be printed 

out and the program will terminate. 
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4 . Warning Concerning Use of Ship Assignment Doctrine 2 



Assignment doctrine 2 has one definite shortcoming. 
It only considers which types of ships are present in the 
problem, not \>fhich types are available to fire at the time 
the mission is called for. Thus, a situation could occur 
in which the only types of ships appropriate to fire a 
mission are off-station for rearming or refueling at the 
time the mission is called for. Then the mission would be 
required to wait until one of these ships returns before it 
is fired although other ships with guns of adequate range 
might be available probably, resulting in an unacceptably 
large amount of queue time. 

E. SHIP AVAILABILITY 

After a mission has been assigned a mission type in ac- 
cordance with the rules of the specified assignment doctrine, 
subroutine AVAIL is called once for each type of ship which 
is eligible to fire the mission and present in the problem. 
The arguments of AVAIL are the number of ships present of 
the type being considered, and the type of ship being con- 
s idered . 

1 . Ship Status ^ 

The status of each ship in the problem is maintained 
by the variable STATUS(S,T), where S is the ship number and 
T is the ship type. There are four possible values of 
STATUS(S,T), as follov,rs: 

1: Ship available for assignment. 

2: Ship presently assigned to a mission. 
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3: Ship presently refueling, or in a "refuel after 

mission" status, or engaged in a combined refueling/rearming. 

4: Ship presently rearming. 

2 . Vector of Available Ships 

The series of calls on AVAIL for each mission re- 
sults in a vector (VAIL) of ships of types appropriate for 
the mission for which STATUS(S,T) = 1. The values in this 
vector are two-digit numbers with the ship number as the 
first digit and ship type as the second digit. For example, 
VAIL(l) = 23 means that the first ship checked which was 
available was 5"/54 ship #2. S/R AVAIL also keeps a tally 
of the number of ships which are available to fire the mis- 
sion. This tally is the value of the variable KOUNT . 

3 . Assignment of Closest Ship 

After the vector of available ships is obtained by 
calls on AVAIL, control is returned to the main program. If 
KOUNT = 0, the mission joins the queue. If KOUNT > 0, then 
the distance of each ship in the vector VAIL from the firing 
position is determined in the following manner: 

Each available ship is considered to be at some 
random point in the ^^raiting pattern. For each ship, its 
distance from the firing position is calculated as a uni- 
formly distributed random variable between the minimum pos- 
sible distance (MINDIS) and the maximum possible distance 
(MAXDIS) . The method of calculation of MINDIS and M\XDIS 
is indicated in Figure 2. 

The closest ship is then assigned to fire the mis- 
sion. This is done by the following calculations: 
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TYPSHP(I) = MOD (VAIL(JD), 10) 
HULLNRCI) = VAIL(JD)/10 



where JD is the position in the vector VAIL of the closest 
ship . 

F. THE QUEUE 

The queue in the model consists of those missions for 
which the time of call of fire is less than the current prob- 
lem time (MSNSKD(I) < CLOCK) but to which a ship has not yet 
been assigned. Thus, missions in the queue are referred to 
as "awaiting assignment". 

1 . Joining the Queue 

A mission joins the queue when no ship appropriate 
for the mission type is available at the time the mission is 
called for'(KOUNT = 0). Mission numbers of missions awaiting 
assignment are maintained in the vector Q. The current num- 
ber of missions awaiting assignment is maintained by the vari- 
able INQ. In view of the assumption stated in paragraph II. 

D. 5 , the maximum possible number of missions in the queue is 
equal to the number of patrol subareas. 

2 . Leaving the Queue 

A mission leaves the queu^ when an appropriate type 
ship becomes available by having completed another mission or 
having completed refueling and/or rearming, and there are no 
other missions ahead of it in the queue which are waiting 
for the same type of ship. Each time a ship becomes avail- 
able and INQ > 0, each mission in the queue, beginning with 
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Q(l) is checked to see if the available ship is appropriate 
to that mission. As soon as one is located, that mission is 
removed from the queue and all other missions following it 
are moved up one number in the queue. 

3 . Queue Discipline 

There are two types of queue discipline available 
to the user. The type desired is specified by the input 
variable QDISC, which may have the following values: 

1: First-come first served. 

2: Personnel missions have priority. 

The use of QDISC = 2 places type 1 and 2 missions ahead of 
type 3 and 4 missions in the queue. There is no priority 
relationship between type 1 and 2 missions. Note that 
QDISC = 2 should not be used when ASMDOC = 2, or vice-versa, 
since ASNDOC = 2 only allov'/s 5”/38 ships to fire personnel 
missions. There is no provision in the program for missions 
which are already in progress (i.e., a ship has been assigned) 
to be preempted by any other mission. Thus, once a ship is 
assigned to fire a mission, it is not available for any other 
mission until that mission is completed. 

G. SERVICE TIME 

Each time a ship is assigned to fire a mission, subrou- 
tine SERVIC is called to calculate the service time for that 
mission. The arguments of SERVIC are type of ship assigned, 
target type, target range, mission number, and distance of 
assigned ship from the firing position. S/R SERVIC calcu- 
lates the total service time of the mission and amount of 
ammunition expended during the mission. 



1 . Total Service Time CSERVC(I)) 

There are three components of total service time: 
a. Travel Time (TRAVEL (I)) 

Travel time is the time required for the ship 
to reach the firing position after it has been assigned. It 
is calculated in one of three ways. 

(1) Ship Coming from, the Waiting Pattern . If 
the assigned ship was assigned by virtue of having been the 
closest available ship of appropriate type at the time the 
mission was called for, then its travel distance is the dis- 
tance calculated by the method described in Section III.E.3. 

(2) Ship Having Just Completed Another Mission . 
In this case the travel distance is the distance along the 
gunline between the two firing positions. This distance is 



RUN = lABS (EWTAR(E) -EWTAR(P)) 



where E is the number of the mission just completed and P is 
the number of the new mission to which the ship has been as- 
signed. 

(3) Ship Having Just Completed Refueling and/or 
Rearming . A ship is considered to have completed refueling 
and/or rearming when it returns to a position midway between 
the northern and southern edges of the waiting pattern and 
due south of the firing position. Then travel distance is 

RUN = WTDIS + NSWAIT/2. 



Regardless of which of the above methods is used 
to calculate travel distance, S/R SERVIC calculates travel 
time as 
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TRAVEL(I) = DIS/CSPEED * 33.33) 



where DIS is the appropriate travel distance, SPEED is an 
input variable which specifies the speed (in knots) at which 
ships proceed to the firing position, and 33.33 is a factor 
which converts knots (nautical railes/hour) to yards/minute. 

b. Setup Time (SETUP(I)) 

This is the time required for a ship to fix its 
navigational position and ready its guns after it has reached 
the firing position. S/R SERVIC calculates setup time as a 
uniformly distributed random variable between the minimum ex- 
pected setup time (MINSET) and maximum expected setup time 
(MAXSET) , where MINSET and MAXSET are input variables spec- 
ified by the user. 

c. Firing Time (FIRE(I)) 

Actual firing time is the time interval between 
the time the first round is fired and the time the last round 
is fired. For each mission this is considered to be an ex- 
ponentially distributed random variable with a mean which is 
a function of the type of ship firing, the target type, and 
the target range. 

(1) Mean Firing Time Due to Target Range . The 

% 

portion of the mean firing time which is due to the target 
range is considered to be a linear function of the target 
range. This portion is calculated’ as RANGE (I) /DENOM (J , K) , 
where J is the type ship and K is the target type. The 
values of DENOM(J,K) are specified by the user for each pos- 
sible combination of ship type vs. target type. Thus, DENOM 
is a 4 X 4 matrix. 
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(2) Mean Firing Time Due to Target Type. This 
quantity is specified by the user as FIRTIM(J,K) where J is 
the ship type and K is the target type. Thus, FIRTIM is 
also a 4 X 4 matrix. 

The mean of the exponentially distributed firing 
time for each mission is calculated as 

MU = RANGE(I)/DENOM(J,K) + FIRTIM(J,K). 

Therefore, the firing time for each mission is calculated as 

FIRE(I) = -MU*ALOG(XRN) . 

After travel, setup, and firing time have been 
calculated, total service time is calculated as 

SERVC(I) = TRAVEL(I) + SETUP(I) + FIRE(I). 

2 . .Ammunition Expenditure (4ALMEXP(I)) 

S/R SERVIC also calculates the amount of ammunition 
expended during each mission. The user specifies the minimum 
(MINAMM(J , K) ) and maximum (MAXAMM(J , K) ) expected rounds of 
ammunition expended per minute by ship type J firing against 
target type K. These values should take into account the 
average rounds per minute over an entire mission, i.e., the 
slow rate of fire at the beginning of the mission and the 
rapid rate of fire after the fire has been spotted onto the 
target should be averaged. Therefore, MAXA^tM(J,K) should 
normally not be the maximum attainable rate of fire by a 
type -J ship . 

The number of rounds of ammunition expended per minute 
for each mission is calculated as a uniformly distributed 
random variable between MINAMM(J,K) and MAXAMM(J,K). Therefore 
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the number of rounds of ammunition expended during the Ith 
mission is calculated as 

AMMEXP(I) = (MINAMM(J,K)+XRN*(MAXAMM(J,K) -MINAT-IM ( J , K) ) ) *FIRE ( I) 
where J is the type of ship and K is the type of target. 

H. REFUELING AND REARMING 

The user has two options available for taking refueling 
and rearming of ships into account in the problem. The op- 
tion desired is specified by the input variable REPPOL (Re- 
plenishment Policy), which has values defined as follows: 

1: Ships do not refuel or rearm. 

2: Ships refuel according to a predetermined schedule 

and rearm when their on-board ammunition level reaches a 
specified minimum level. 

The case in which ships do not rearm or refuel would 
arise when enough ships are available to replace ships which 
are required to depart for refueling and/or rearming, thus 
the departure of the affected ship would not effect the 
queueing results. 

If REPPOL = 2, refueling and rearming is handled in the 
following ways in the program: 

1 . Refueling 

a. Methods of Establishing a Refueling Schedule 
The user establishes a refueling schedule for 
all ships for the entire duration of the problem. There are 
two options available for establishing this schedule. The 
desired option is specified by the input variable REFSKD, 
which has values defined as follows: 



1: The user specifies the exact time of each 

scheduled refueling for each ship. For each type of ship, a 
data card is prepared which contains the floating point time 
(in days from zero problem time to 2 decimal places] that 
each ship of that type is to depart station for its first 
scheduled refueling. Similar cards are prepared for the 
2nd, 3rd, and up to the 9th scheduled refuelings for ships 
of each type. These refueling times are the values of the 
variable XFULSK(I , J,K] , where I is the ship number, J is 
the ship type, and K is the number of the refueling. The 
total number of rounds of refueling is specified by the in- 
put variable FINFUL (must be 59). Thus, there must be 4 x 
FINFUL cards in the refueling schedule data deck . The pro- 
gram converts these scheduled refueling times to minutes 
from zero problem time and places them on the Event Chain in 
the proper chronological order. 

2; The refueling schedule scheme discussed 
above is replaced by a simpler (although less flexible) 
scheme. The entire schedule is specified on one data card. 
The first five quantities on the card are the floating point 
times (in days from zero problem time to 2 decimal places) 
that the first cruiser is to deparf for its 1st thru 5th 
scheduled refueling. These are specified by the input vari- 
ables XCRFST(I) , 1=1 thru 5. The next 5 quantities are the 
times the first destroyer is to depart for its 1st through 
5th scheduled refuelings. These are specified by the input 
variables XDDFST(I) , 1=1 thru 5. The last two numbers are 
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the time interval in minutes at which succeeding cruisers 
and destroyers follow the first for refueling. These are 
specified by the input variables CRINT and DDINT. The pro- 
gram then calculates a refueling schedule for the entire 
duration of the problem using these parameters. There are 
several rules which must be followed in the preparation of 
the refueling schedule data cards. 

a. For ships which are not present in the 
problem, the appropriate columns on the data cards should be 
left blank. 

b. If REFSKD = 1, for ships which are present 
in the problem but are not scheduled to refuel FINFUL times, 
a large number (99.00 for example) must be inserted in the 
appropriate columns on the data cards. 

c. If REFSKD = 2 and less than five rounds of 
refueling are to be used, a large number (99.00) must be in- 
serted in the appropriate columns on the data card. 

d. If REFSKD = 1 and no ships of a particular 
type are present in the problem, a blank card must be in- 
serted in the data deck in the appropriate place for each 
round of refuelings. 

Examples 

REFSKD = 1 : The expected duration of the problem is 10.0 

days. NR855 = 2, NR647 = 0, NR554 = 1, NR538 =2. It is 
desired that the cruisers refuel on the sixth day and the 
destroyers refuel on the fifth and ninth days. Then typical 
values of XFULSK(I , J , K) would be 
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First Card 



Fifth Card 



XFULSK(1 ,1 ,1) = 5.50 


XFULSKCl ,1 ,2) = 99.00 


XFULSK(2 ,1,1) = 5.75 


XFULSKC2,1,2) = 99.00 


Third Card 


Seventh Card 


XFULSKCl ,3,1) = 4.25 


XFULSKCl ,3,2) = 8.25 


Fourth Card 


Eighth Card 


XFULSKCl, 4,1) = 4.50 


XFULSKCl, 4, 2) = 8.50 


XFULSK(2,4,1) = 4.75 


XFULSKC2,4,2) = 8.75 



FINFUL = 2 (2 rounds of refueling) . Therefore there would - 
be 8 data cards, of which the 2nd and 6th data cards would 
be blank since NR647 = 0. 

REFSKD = 2 : For the same situation described above, the fol- 

lowing values on one data card would produce the same schedul 
as above : 

XCRFSTCl) = 5.50 XDDFST(l) = 4.25 

XCRFST(2) thru XCRFST(5) = 99.00 XDDFST(2) = 8.25 

CRINT = 360 XDDFST(3) thru 

XDDFST(5) = 99.00 

DDINT = 360 

b. Latest Fueling Time 

The program calculates the time of extraction 
of each patrol from its subarea as INSERT(NR) plus a uni- 
formly distributed random variable between XSHORT(NR) and 
XLONG(NR) . The latest of these extraction times is stored 
as the variable MAXTRC. In computing the fueling schedule 
according to user input parameters, no refuelings are placed 
on the Event Chain which have a scheduled time later than 
MAXREF = MAXTRC - NOMO , where NOMO is an input variable 



specified by the user. MAXREF is defined as the latest time 
that ships will depart for refueling. For example, if it 
is desired that no ship depart for refueling later than 6 
hours before the latest anticipated extraction time, then 
NOMO = 360. 

2 . Rearming 

The user specifies several input variables which 
are used to determine when each ship requires rearming. These 
input variables are defined as follows: 

INAMM(S,T): The initial ammunition level (rounds 

on board) of ship number S of type T. 

MINA>LMCJ,K) and MAXAMM(J,K) : Defined in Section 

III.G.2. 

MINLEV(T): The minimum ammunition level to v^hich a 

ship of type T is allo\\red to fall before it is required to 
rearm. 

Each time a mission is fired, ammunition expenditure 
for that mission is calculated by S/R SERVIC as discussed in 
Section III.G.2. This amount is subtracted from the current 
ammunition level of the ship which fired the mission, as fol- 
lows : 

AMMLEV(S,T) = AMMLEv\s,T) - AMMEXPCI). 

If AMMLEV(S,T) < MINLEV(T), the ship is sent to rearm. 

3 . Rearming and Refueling Completions 

a. Refueling 

When an event for which EVTYP(K) = 3 (scheduled 
refueling) occurs, the status of the scheduled ship is checked. 



I£ the ship is presently firing a mission (STATUS(S,T) = 2) 
it is placed in a "refuel after mission" status (STATUS(S,T) 
=3). If STATUS (S,T) = 1, the ship is sent to refuel im- 
mediately. In either case, the time of completion of fueling 
is calculated as follows: 

The user specifies the minimum and maximum ex- 
pected lengths of time that it should take for a cruiser and 
a destroyer to refuel. These are the input variables MINFLC, 
MAXFLC, MINFLD, and MAXFLD. These times should take into 
account travel time to and from the rendezvous point as well 
as time alongside the oiler. For example, if it is planned 
that ships will rendezvous with the oiler at a point 20 
miles from the center of the waiting pattern and will proceed 
to the rendezvous point at a speed of 20 knots, and that a 
destroyer should require 30 to 60 minutes alongside the 
oiler, than MINFLD should be 60+30+60 = 150, and MAXFLD 
should be 60+60+60 = 180. 

The time of completion of refueling (COMFUL) is 
calculated as CLOCK time (the time the ship departed for re- 
fueling, either after a mission is completed or at the sched- 
uled refueling time) plus a uniformly distributed random 
variable between MINFLC and MAXFLC for cruisers, or MINFLD 
and MAXFLD for destroyers. COMFUL is then placed on the 
Event Chain. 

b. Rearming 

Completion of rearming (COMARM) is handled in 
the same way as COMFUL except that the parameters of the 
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uniform distribution of rearming time are MINARC, MAXARC, 
MINARD, and MAXARD, specified by the user with the same 
considerations as the refueling time parameters. When 
COMARM is calculated, the rearming ship's ammunition level 
(AMMLEV(S ,T) ) is set equal to its initial ammunition level 
(INAMM(S,T)) . 

c. Combined Refueling/Rearming 

When a ship which is scheduled to refuel is 
off-station for rearming at the time the scheduled refueling 
occurs, then a refueling time calculated in the manner de- 
scribed above is added to the previously calculated COMARM. 
From this time, the input variable COMBO is subtracted. 

COMBO is defined as the amount of time off-station saved by 
conducting both refueling and rearming in the same evolu- 
tion. Thus, it should be approximately equal to the time 
for a one-way trip from the ammunition ship to the waiting 
pattern plus the time for a one way trip from the waiting 
pattern to the oiler, minus the time required to travel from 
the ammunition ship to the oiler. COMBO must not be greater 
than MINFLD or MINFLC . 

The previously calculated C0MAR?4 is removed 
from the Event Chain and the new time of completion of com- 
bined refueling/rearming (COMREP) is placed on the Event 
Chain . 

I . OUTPUT DATA 

All output of the program is well - formatted to facili- 
tate reading and interpretation. Each type of output is dis- 
cussed belcw. 



1 . Printout of Input Data 



A printout of the user-specified values of most of 
the input variables is produced at the beginning of the out- 
put data. Appendix C contains an example of this printout. 

2 . Mission History 

The Mission History is a chronological list of the 
12 most significant parameters assigned to each mission dur- 
ing a particular iteration of the problem. There are three 
options available to the user concerning the Mission History 
printout. The desired option is specified by the input 
variable MISHIS, which has values as follows: 

1: No Mission History printout desired. 

2: Printout Mission History for the first iteration 

only . 

3: Printout Mission History for each iteration. 

The time of call for fire and completion of each mission are 
converted to date-time groups in the Mission History printout. 
The first 1 or 2 digits are the problem day and the last four 
digits are the 24 -hour clock time. 

Appendix D contains an example of the Mission His- 
tory printout. 

3 . Statistical Results for Edch Iteration 

For each iteration, a printout of the pertinent 
queueing, ship utilization, refueling, rearming, and ammu- 
nition expenditure statistics is available. These are the 
results for that iteration only, they are not cumulative 
over all iterations up to that point. The user specifies 



58 



whether printout of iteration statistical results is desired 
by the input variable ITSTAT, which may have the following 
values : 

1* Printout of iteration statistical results is 

desired . 

2: Printout of statistical results for each itera- 

tion is not desired. 

Appendix E contains an example of the printout of 
iteration statistical results. 

4 . Overall Statistical Results 

After all NRIT iterations of the problem have been 
completed, a printout of the pertinent overall statistical 
results is produced. This printout includes the mean and 
the standard deviation of the mean for each quantity except 
for refueling, rearming, and ammunition expenditure data, 
for which the mean only is printed out. 

In computing the mean and the standard deviation 
of the mean of the quantities in the overall statistical re- 
sults, the following formulas are used: 




where x^ is the result of the i^h iteration and n = NRIT. 

Using the means and standard deviation from the 
overall statistical results for NRIT iterations, confidence 
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intervals for the actual value of each output quantity can 
be calculated using the following formula: 



X 






t(n-l) 

(l-a)/2 



< U < X + 0= 



tCn-1) 

Cl-a)/2 



where x is the observed mean, o= is the observed standard de- 

viation of the mean, and t^l? is the value from a table 

(l-a)/2 

of values of the Student’s-t distribution with n-1 degrees 
of freedom at the (l-a)/2 level of significance. a is de- 
fined as the probability that the actual value of the quan- 
tity lies outside of a confidence interval calculated by 
the above formula. Therefore (1-a) is the degree of con- 
fidence with which statements concerning the actual value 
of the quantity can be made. 

Example : 

U = AVERAGE WAITING TIME PER MISSION 

10 iterations of a problem yield x = 20.0 and a= = 1.00. The 
user desires to be able to predict with 95% confidence the 
interval in which the actual value of AVERAGE WAITING TIME 
PER MISSION lies based on these results of 10 iterations. 

A table of values of the Student 's-t distribution shows 
t ®972 = 2.26. Therefore, it can be stated with 95% con- 
fidence that 



20.0 -(1.00) (2.26) < y < 20.0 + (1 . 00 ) (2 . 26 ) 
or 

17.74 < AVERAGE WAITING TIME PER MISSION < 22.26. 
Appendix F contains an example of the printout of 
overall statistical results. 
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J. INPUT DATA - DATA DECK 

This section lists all input variables in the data deck 
by card and column number(s) . For each variable the follow- 
ing information is listed: variable name, subsection of 

section III in which the variable is described, the FORTRAN 
variable type, the maximum possible value which can be as- 
signed, and the units of the indicated quantity. All quanti 
ties must be right- adj usted within their appropriate columns 

1. First Card: Miscellaneous Problem Definition Vari- 



Cols 


ables 

Variable 


Re fe rence 


IZ£g- 


Maximum 


Units 


1-2 


NRAREA 


B.l 


12 


50 


Sub area; 


3 


NR855 


A. 5 


11 


9 


Ships 


4 


NR64 7 


A. 5 


11 


9 


Ships 


5 


NR554 


A. 5 


11 


9 


Ships 


6 


NR5 38 


A. 5 


11 


9 


Ships 


7-11 


MAX855 


D. 3 


15 


99999 


Yards 


12-16 


MAX647 


D. 3 


15 ■ 


99999 


Yards 


17-21 


MAX554 


D. 3 


15 


99999 


Yards 


22-26 


MAX538 


D. 3 


15 


99999 


Yards 


27-31 


MAXNGF 


D. 3 


15 


99999 


Yards 


32-38 


WTNW 


B.3 


17 


9999999 


Yards 


39-45 


WTNE 


B.3 


17 ‘ 


9999999 


Yards 


46-51 


NS WAIT 


B.3 


16 


999999 


Yards 


52-57 


WTDIS 


B.3 


16 


999999 


Yards 


58-59 


SPEED 


G. 1 . a 


12 


99 


Knots 


60-61 


MINSET 


G.l.b 


12 


99 


Minutes 


62t63 


MAXSET 


G. 1 . b 


12 


99 


Minutes 


64-67 


HHOUR 


C.2 .a. (1) 


14 


1439 


Minutes 



2. 


Second Card 


. : Program 


Control 


Indicators 




Cols 


Variable 


Re Terence 


Type 


Maximum 


Units 


1 


ASNDOC 


D.l 


11 


2 


— 


2 


QDISC 


F.3 


11 


2 


— 


3 


REPPOL 


H 


11 


2 


— 


4 


REFSKD 


H.l 


11 


2 


— 


5 


FINFUL 


H.l 


11 


9 


Rounds of 












Re fueling 


6 


MI SHIS 


1.2 


11 


3 


— 


7 


ITSTAT 


1.3 


11 


2 


— 


8-10 


NRIT 


A. 7 


13 


999 


Iterations 


3. 


Next NRARE;\ 


. Cards: Subarea Information 


Cards 




One card is 


prepared 


for each 


patrol subarea. Cards 


are placed in the data deck in 


order of subarea n 


lumbers. NR 


is the 


subarea number. 








Cols 


Variable 


Reference 


Type 


Maximum 


Units 


1-7 


NORTH (NR) 


B.2 


17 


9999999 


Yards 


8-14 


SOUTH (NR) 


B.2 


17 


9999999 


Yards 


15-21 


EAST (NR) 


B.l 


17 


9999999 


Yards 


22-28 


WEST (NR) 


B.l 


17 


9999999 


Yards 


29-33 


XCFF(NR) 


C . 2 . a 


F5 .1 


999.9 


Missions 


34-38 


XSHORT(NR) 


C . 2 . a 


F5.1 


999.9 


Days 


39-43 


XLONG(NR) 


C . 2 . a 


F5. 1 


999.9 


Days 


44-49 


INSERT (NR) 


C.2 .a. (1) 


16 


999999 


Minutes 


50-54 


XPROBl(NR) 


C.2 .b. (5) 


F5 . 2 


1.00 


Probab ility 


55-59 


XPROB2 (NR) 


C.2.b. (5) 


F5.2 


1.00 


Probability 


60-64 


XPROB3(NR) 


C.2.b. (5) 


F5.2 


1.00 


Probability 


65-69 


XPROB4(NR) 


C.2.b. (5) 


F5 . 2 


1.00 


Probability 
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4. Next 2 Cards: Parameters of Mean Firing Time 



First Card: Mean Firing Time Due to Target Type 



Cols 


Variable Reference 


Type 


Maximum 


1-2 


FIRTIM(1 ,1) G.l.c. (2) 


12 


99 


3-4 


FIRTIM(2,1) " 


ft 


ft 


5-6 


FIRTIM(3,1) " 


ft 


f ? 


7-8 


FIRTIM(4,1) 


ft 


ft 


9-10 


FIRTIM(1,2) 


ft 


f f 


11-12 


FIRTIM(2,2) 


ft 


ft 


13-14 


FIRTIM(3,2) 


ft 


ft 


15-16 


FIRTIM(4,2) 


ft 


ft 


17-18 


FIRTIM(1,3) 


ft 


ft 


19-20 


FIRTIM(2,3) 


ft 


ft 


21-22 


FIRTIM(3,3) 


ft 


ft 


23-24 


FIRTIM(4,3) 


ft 


ft 


25-26 


FIRTIM(1,4) 


ft 


ft 


27-28 


FIRTIM(2,4) 


ft 


ft 


29-30 


FIRTIM(3,4) 


ft 


ft 


31-32 


FIRTIM(4,4) 


ft 


ft 


Second 


Card: Denominator of 


Range Component 


1-5 


DENOM(l,l) G.l.c. (1) 


15 


99999 


6-10 


DENOM(2,l) 


ft 


ft 




• 

(Same order of subscripts as 


FIRTIM) 


76-80 


• 

DENOM(4,4) 


ft 


ft 



(IF REPPOL = 1, skip to section III.J.7) 



Units 

Minutes 

TT 

M 

ft 

tf 

TT 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

Yards/Minute 

ft 
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Refueling Data 



Cols 

1-4 

5-8 



a. First Card: Miscellaneous 

Variable Reference Type 

NOMO H.l.b 14 

COMBO H . 3 . c 14 



Maximum Units 
9999 Minutes 
9999 Minutes 



(IF REFSKD = 2, skip to section III.J.S.c) 



b. Next 4 X FINFUL Cards: User- Input Refueling 

Schedule ~ 



First Card 



1-7 


XFULSK(1 ,1 ,1) 


H.l 


VI .2 


9999 


8-14 


XFULSK(2,1 ,1) 


tt 


f 1 


tt 


57-63 


XFULSK(9,1,1) 


ft 


ft 


tt 


Second 


Card 








1-7 


XFULSK(1,2,1) 


t1 


tt 


:t 


8-14 


XFULSK(2,2 ,1) 


M 


tt 


tt 


57-63 


XFULSK(9,2,1) 


M 


tt 


tt 


Third ' 


Card 








1-7 


XFULSK(1 ,3,1) 


If 


tt 


tt 


8-14 


XFULSK(2 ,3,1) 
• 


ft 


tt 

\ 


tt 


57-63 


XFULSK(9,3,1) 


tf 


tt 


tt 


Fourth 


Card 









1-7 XFULSK(1 ,4 ,1) 

8-14 XFULSK(2 ,4,1) 



XFULSK(9,4,1) 



57-63 



' ® 





Next 4 Cards 



Same as first four cards except that third subscript is 2, 



Continue in the same manner until there are FINFUL sets of 



4 cards each. There 


must be 4 X 


FINFUL 


cards in 


the user- 


input 


refueling schedule section 


of the 


data deck. 


(Skip 


to section III 


. J. 5 .d) 










c. Next 


Card: Parameters of 


Computer- 


Calculated 


Refueling Schedule 












Cols 


Variable 




Reference 


Type 


Maximum 


Units 


1-6 


XCRFST(l) 




H.l .a 


F6.2 


999.99 


Days 


7-12 


XCRFST(2) 




ft 


It 


tt 


1 1 


13-18 


XCRFST(3) 




It 


tt 


tt 


1 1 


19-24 


XCRFST(4) 




»l 


It 


tt 


1 1 


25-30 


XCRFST(5) 




M 


It 


1 1 


tt 


31-36 


XDDFST(l) 




If 


It 


tt 


tt 


37-42 


XDDFST(2) 




M 


tt 


tt 


It 


43-48 


XDDFST(3) 




tl 


tl 


1 1 


tt 


49-54 


XDDFST(4) 




It 


tl 


tt 


tt 


55-60 


XDDFST(5) 




tl 


tt 


1 1 


tt 


61-65 


CRINT 




1 1 


15 


99999 


Minutes 


66-70 


DDINT 




tl 


15 


99999 


Minutes 




d. Next 


Card : Minimum 


and Maximum Refueling Off 




Station 


Times 








1-5 


MINFLC 




H.3.a 


15 


99999 


Minutes 


6-10 


MAXFLC 




It 


tt 


. tt 


1 1 


11-15 


MINFLD 




1 1 


tt 


Tt 


tt 


16-20 


MAXFLD 




It 


It 


tt 


It 







i 
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. Rearming Data 

a . First Card: Minimum Allo\ved Ammunition Levels 

and Minimum and Maximum Rearming Off-Station 



Cols 


Times 

Variable 


Reference 


Type 


Maximum 


Units 


1-6 


MINLEV(l) 


H.2 


16 


999999 


Rounds 


7-12 


MINLEVC2) 


ft 


It 


1 1 


II 


13-18 


MINLEV(3) 


tl 


It 


II 


II 


19-24 


MINLEV(4) 


1 1 


1 1 


It 


II 


25-29 


M INARC 


H.3.b 


15 


99999 


Minutes 


30-34 


MAXARC 


tl 


II 


It 


II 


35-39 


MINARD 


tl 


II 


1 1 


1 1 


40-44 


MAXARD 


II 


II 


1 1 


II 



b. Next 4 Cards: Initial Ammunition Levels 



First Card 



1-7 


INAMMCl,!) H.2 


17 


9999999 


Rounds 


8-14 


INAMM(2,1) 


It 


1 1 


1 1 



57-63 INAMM(9,1) 

Second Card 
1-7 INAMM(1,2) 

8-14 INAMM(2,2) 

• • 

57-63 INAMM(9,2) 

Third Card 

1-7 INA:vIM(1,3) 

8-14 INA^LM(2,3) 

57-63 INAMM(9,3) 



Fourth Card 



1-7 


INAMM(1,4) H.2 


17 


9999999 


Rounds 


8-14 


INAMM(2,4) 

• 


ft 


tr 


ft 


57-63 


• 

INAiMM(9,4) 


f f 


ft 


ft 



7 . Next 2 Cards : Ammunition Expenditures Per Minute 



Cols 


Variable 


Reference Type Maximum Units 


First 


Card 




1-4 


MINAMM(1 ,1) 


G.2 14 9999 Rounds/Minute 


5-8 


MINAMM(2 ,1) 


ff If ff ff 




(Same order 


of subscripts as FIRTIM) 


61-64 


MINAiMM(4 ,4) 


M If ff ff 


Second 


Card 




1-4 


MAXME.1(1 ,1) 


ff ff ff ff 


00 

1 • 

LO 


MAXAiMM(2 ,1) 


ff ff ff ff 




(Same order 


of subscripts as FIRTIM) 


61-64 


MAXAMM(4 ,4) 


ff ff ff ff 


8 


. Last Card; 


Random Number ^Generator Initializers 




This card contains a string of 70 random digits in 



columns 1-70, with the requirement that every fifth column 
contains an odd number. 



K. COMPUTER STORAGE AND TIME REQUIREMENTS 

The program runs on an IBM- 360/67 at the W. R. Church 
Computer Center at the U. S. Naval Postgraduate School. For 
other systems, the storage and time requirements may differ 
from those stated below. 

1 . Computer Storage 

The amount of core area required is dependent on the 
amount of array storage set aside for array variables, which 
of course depends on the estimated size of the problem to be 
simulated. The program as presented in Appendix B allows up 
to 2300 events and 1000 missions to be generated during each 
iteration of the problem. Thus the 19 parameters associated 
with each mission all have dimension size 1000 and the 4 pa- 
rameters associated with each event all have dimension size 
2300. With these array sizes, the program requires approxi- 
mately 195K bytes of execution core area. 

2 . Computer Central Processing Unit Time 

The program requires approximately 1 minute and 25 
seconds of CPU time for the compile and link steps. CPU time 
required in execution is primarily dependent on the number 
of missions generated and the number of iterations performed. 
The amount of execution time can be calculated approximately 
as 1/200 seconds per mission per iteration. The expected 
number of missions per iteration is approximately^ XCFF (NR) . 

Therefore, if the expected number of missions per 
iteration for a problem is 500 and 10 iterations are perform- 
ed (NRIT = 10) , then 



Execution Time = 500 (i/200) (10) = 25 seconds. 



IV. VERIFICATION AND VALIDATION OF THE MODEL 



A. VERIFICATION 

The fact that logic flow in the program is as intended 
and that output statistics accurately describe the situa- 
tion being simulated has been verified in several ways. 

These are discussed below. 

1 . Verification of Intermediate Calculations and Data 

Flow 

During the formulation of the program and numerous 
debugging runs, results of calculations and other pertinent 
data at various points in the program were printed out to 
determine their correctness. When the program was considered 
finally debugged, such checks were performed at all points 
in the program at which calculations were performed or sig- 
nificant data points referenced. The results of these 
checks were carefully examined, and all values were deter- 
mined to be correct. 

2 . Tracing 

Logic flow through all possible branches of the 
program has been carefully traced and determined to be as 
intended . 

3 . Extreme Value Problems 

When the program was considered finally debugged, 
several problems were formulated using extreme foreseeable 
maximum and minimum values of input parameters. In each 
case, the program successfully ran to completion with 
statistical results reflecting the extreme values. 



69 




1 




4 . Verification of Exponential Distribution 



To verify that mission interarrival times are 
actually produced according to the exponential distribution, 
a Kolmogorov-Smirnoff test for goodness-of-fit was perform- 
ed on the interarrival times in a subarea picked at random 
in a trial run of the program. The subarea picked was sub- 
area 11. The following input data was used for subarea 11 
in the trial run: 



XCFF(ll) = 35.0 missions 
XSHORT(ll) = 6.5 days 

XLONG(ll) = 7.5 days 

Thus , 



XDURAT(ll) = 7.0 days 

XCALRTCll) = .00347 missions/minute. 

The continuous curve in Figure 3 is a plot of the 
CDF of an exponential distribution with rate = .00347. The 
step function in Figure 3 is the sample CDF of the 23 mis- 
sion interarrival times generated for subarea 11 in the 
trial run. 

The Kolmogorov-Smirnoff test [pp . 238-240, Ref. 5] 
consists of observing the absolute difference between the 

k 

value of the sample CDF and the value of the hypothesized 
CDF at all points in the sample. The test statistic is 



D^ = max { I F(x) - F^^Cx) | } 
n 

where F(x) is the value of the hypothesized CDF at x and 
Fj^(x) is the observed value of the sample CDF at x for a 
sample of size n. 



f 



i 




Figure 3. Kolmogorov -Smirnoff Test for Goodness -o f- Fit 

of Mission Interarrival Times to the Exponential 
Distribution . 
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Inspection o£ Figure 3 yields ^2'S ~ Table 

VI in Reference 5 shows that < .276 is acceptable with 

95% confidence. Therefore the sample of interarrival times 
in subarea 11 successfully passes the' Kolmogoro v-Smirnof f 
test for goodness -of - fit for the exponential distribution. 
From this result, it was assumed that all interarrival times 
and mission firing times fit an exponential distribution, 
since the same method of generating exponentially distributed 
random variables was used in each case. 

5 . Comparison with GPSS Model 

A short program in the General Purpose Simulation 
System (GPSS) was written to simulate the system described 
herein in a very basic mode, which was easily adaptable to 
GPSS. The program described herein was run using input 
values which placed it in exactly the same mode as the GPSS 
program. The queueing results of each were then compared 
and found to be in very close agreement. 

6 . Verification of Relationships Between Output 

Statistics 

By investigating the mathematical relationships be- 
tween various output statistics, it was verified that the 
calculations which produce the output values are properly 
performed. An example of an investigation of such a rela- 
tionship is discussed belo\\?. 

OVERALL SHIP UTILIZATION: If refueling and rearming 

are not used in the problem, overall ship utilization should 
equal 
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Total Missions x Average Service Time Per Mission 
Number of Ships x Total Problem Time 

For a trial run of the program in which refueling and rearm- 
ing were not used, the folloi\ring results were obtained: 

NR855 = NR647 = NR538 = 1 
NR554 = 0 

TOTAL NUMBER OF MISSIONS FIRED = 385 
AVERAGE SERVICE TIME PER MISSION = 43.0 minutes 
TOTAL PROBLEM TIME = 8 days , 15 hours , 54 minutes 
OVERALL SHIP UTILIZATION = 0.44. 

Using these figures in the formula stated above yields: 

Overall Ship ^ 385 x 4 3.0 _ - . 

Utilization 3 x (8 x 1440 + 15 x 60 + 54) 

Since overall ship utilization is not calculated 
using the above formula in the program, but rather by accumu- 
lating an overall sum of minutes of utilization for each 
ship, this result lends further verification to both logic 
flow and accuracy of all of the above listed output statis- 
tics. By methods similar to the above, calculations of 
other related output statistics were also verified. 

7 . Sensitivity Testing 

i ■ 

As a final method of verification, the sensitivity 
of the model to changes in input parameters was tested. An 
example of such tests is discussed below. 

WTDIS: A trial run of the program was made using 

WTDIS = 4000. The output statistics for this problem yielded 
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AVERAGE TRAVEL TIME = 16.5 minutes 



AVERAGE WAITING TIME (ALL MISSIONS) = 8.3 minutes 
PERCENT OF MISSIONS REQUIRED TO WAIT = 26.5%. 

The same problem was then run using WTDIS = 8000, 



all other input values remaining unchanged. It was ex- 
pected that each of these output statistics would increase 
by amounts estimated as follows; 

a. Average travel time should increase by ap- 
proximately 



since each ship coming from the waiting pattern would have 
an additional 4000 yards to travel at a speed of 20 knots, 



coming from the waiting pattern. 

b. Average waiting time (all missions) and percent 
of missions required to wait should increase very slightly 
due to the small increase in average travel time. 

The results of the trial run using WTDIS = 8000 
were as follows : 

AVERAGE TRAVEL TIME = 20.4 minutes 

AVERAGE WAITING TIME (ALL MISSIONS) = 8.5 minutes 

PERCENT OF MISSIONS REQUIRED TO WAIT = 27.2 

These results are very close to those expected, 
further verifying logic flow and program results. 
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B. VALIDATION 



The simulation cannot be considered validated since it 
is obviously infeasible to actually physically conduct a 
problem for experimental purposes within the scope of this 
thesis, and no experimental data of similar exercises is 
available. The design of the model is primarily based on 
personal experience in NGFS both from aboard ship and as a 
spotter ashore, and the results obtained by the simulation 
model appear to be valid within the limits of that experience. 
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V. PARTIAL STDF FEASIBILITY STUDY 



With the computer simulation model finalized and fully 
tested, a partial study of the feasibility of the STDF con- 
cept was undertaken. This study focused only upon the que- 
ueing aspects of providing NGFS from the STDF to SIAF's 
ashore . 

A. OBJECTIVE OF THE STUDY 

The overall objective of the study can be generally 
stated as follows: 

"Determine the optimum number and mix of NGFS ships 
required to adequately support SIAF's in an insur- 
gent-contested area." 

Before proceeding further, the term "adequate" support 
must be defined. Since the model does not consider results 
of missions fired, adequate support cannot be defined in 
terms of results achieved. Rather, the definition must be 
in terms of the queueing aspects of NGFS. 

1 . Decision Criteria 

Two arbitrary decision criteria were established 
in order to define "adequate" suppc^rt. These criteria are 
based solely on personal experience and cannot be considered 
authoritative. The decision criteria are: 

a. The average waiting time per mission should be 
less than 10 minutes. Spotting teams ashore become quite 
anxious and their position becomes increasingly tenuous as 
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the length of time between their call for fire and notifica- 
tion that a ship is on the way increases. Ten minutes seems 
like a reasonable upper bound on the length of time that a 
patrol should have to expect to wait. 

b. Overall ship utilization should be at least 
50%. Paucity of ships to assign to a particular operation 
dictates a requirement such as this, and is the basis of the 
need for a feasibility study of the queueing aspects of NGFS 
of SIAF's. 

2 . Objective Function 

Having established the two decision criteria to 
define "adequate” support, the objective of the study be- 
comes : 

"Determine the minimum number (and mix) of NGFS 
ships required to satisfy the decision criteria." 

In terms of the variables and the output statistics 
in the computer model, the objective can be stated as: 

Minimize: NR855 + NR647 + NR554 + NR538 

Subject to: AVERAGE WAITING TIME (ALL MISSIONS) <10.0 

OVERALL SHIP UTILIZATION > 0.50. 

3 . Optimization Technique 

The optimization technique used can be described 
generally as a 3-step procedure. These steps are: 

a. Hypothesize a typical STDF/SIAF situation. 

b. Exercise the simulation model, varying the 
number of ships of each type which are present. 
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c. Inspect the results of the simulation exercise 
to determine the minimum number (and mix) of ships which 
best satisfies the decision criteria. 

B. HYPOTHESIZED SITUATION - SCENARIO 

In order to exercise the computer simulation, a scenario 
was developed which seems to describe a typical situation as 
envisioned under the STDF/SIAF concept. 

1 . General Situation 

A small country, allied with the U.S., is under 
imminent threat of insurgent takeover. For the past several 
months, there have been numerous reports of infiltration into 
the country and tensions have increased with the country’s 
unfriendly neighbor to the southwest. During the past two 
weeks there have been terrorist attacks in the small fishing 
villages located along the banks of a river which divides 
the country roughly in half. The capitol city is located at 
the mouth of this river. The government now believes that 
the infiltrators are an insurgent force from the unfriendly 
neighboring country, and that the insurgents are building up 
for an attack on the capitol city in an attempt to unseat 
the present government. The allied country has requested 
emergency military assistance from the U.S. to suppress the 
insurgent force. 

2 . Enemy Situation 

The recent terrorist attacks and other reports in- 
dicate that the bulk of the enemy force is dispersed along 
the banks of the river in small guerilla groups. However, 
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since there have been no sightings or use of heavy artillery 
or other significantly large military equipment near the 
river, it is believed that some of the enemy force is en- 
gaged in transporting heavy equipment to strategic locations 
to the west of the capitol city. The most likely routes for 
movement of such equipment are through the sparsely populated 
mountainous regions in the southwest and northwest parts of 
the country. There are estimated to be between 7500 and 
10000 insurgents within the country at the present time. 

3 . Friendly Situation 

The U.S. has agreed to send elements of the STDF 
to assist in suppressing the insurgent force. However, the 
U.S. has made it clear that this military assistance is on 
an emergency basis only, and will be provided for a period 
of seven days only. The STDF commander has received specific 
orders to place only a minimum number of men ashore, and that 
patrols placed ashore are not to actually engage in combat 
except in self-defense. 

a. NGFS/SIAF Detachment 

The STDF commander detached the follo\ving ships 
from the STDF to provide the emergency assistance requested: 

1 LHA (with one Marine Expeditionary Unit 
embarked) 

1 AOE 

X Cruisers and Destroyers 

b. Area of Operations 

(1) Objective Area . The objective area is de- 
fined as the area surrounding the capitol city, 10 miles north 
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and south of the river, and including up to the peaks of the 
mountain ranges in the northwest and southwest. The objective 
area is approximately 20 miles wide and 10 miles deep (inland 
from the coast) . 

(2) Patrol Subareas . Based on expected densi- 
ties of enemy targets and terrain features, the objective area 
has been divided into 34 SIAF patrol subareas. Patrols will 
be inserted into these subareas by helicopter beginning at 
0500 (H-hour) on the first day. Order and time of patrol in- 
sertion will be randomly picked. All patrols will have been 
inserted no later than 2100 on the first day. 

(3) Waiting Pattern . When not assigned to fire 
a mission, ships will operate within a 10 mile wide area 4-6 
miles offshore and centered on the center of the objective 
area since the heaviest target densities are expected to be 
near the center of the objective area. 

c. Firing Doctrine 

In this situation, the U.S. is very much con- 
cerned that its troops ashore not engage in combat unless 
absolutely necessary. Thus it is desired that all calls for 
fire from SIAF patrols be answered as rapidly as possible by 
any ship available, and that personnel targets be fired at 
before material targets if both occur simultaneously. 

d. Refueling and Rearming 

(1) Location of APE . The assigned AOE will re- 
main approximately 15 miles offshore to accomplish refueling 
and rearming. Thus, a one-way trip from the waiting pattern 
to the AOE at a speed of 20 knots should require approximately 
30 minutes. 



80 



(2) Refueling Schedules . Since it is antici- 
pated that the operation will last approximately one week, no 
ship should require refueling more than once. Destroyers are 
scheduled to refuel beginning at 1100 on the fifth day, with 
other destroyers following the first at four-hour intervals. 
Destroyers may also refuel beginning at 1100 on the ninth day 
if the operation lasts that long. Cruisers are scheduled to 
refuel beginning at 1100 on the sixth day, with other cruisers 
following the first at four-hour intervals. It is anticipated 
that destroyers should require 60 to 90 minutes alongside the 
AOE during each refueling, while cruisers should require 90 

to 120 minutes. As the operation progresses and the time of 
withdrawal of all patrols from the objective area becomes 
reasonably firm, no ship will depart station for refueling 
within 12 hours of this anticipated time of total withdrawal. 

(3) Initial and Minimum Ammo Levels . Prior to 
arrival at the area of operations, all NGFS ships have topped 
off their onboard ammunition levels. The STDF commander has 
prescribed minimum allowable ammunition levels for each type 
of ship. These ammunition levels are as follows: 

Type Ship Initial Minimum 



8"/55 

6”/47 

5'754 

5"/38 



8000 

10000 

12000 

6000 



1000 

1000 

2000 

1000 



It is anticipated that destroyers should require 2 to 3 hours 
alongside the AOE to rearm, while cruisers should require 3 
to 4 hours . 
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c. 



HYPOTHESIZED SITUATION - INPUT VARIABLES 



Having hypothesized what appears to be a typical STDF/ 
SIAF situation, the next step was to translate the scenario 
into values of input variables in order to exercise the com- 
puter simulation. Appendix C is the printout of input data 
used in the hypothesized problem, for a run in which NR855 = 
NR647 = NR554 = 1 and NR538 = 0. 

1 . Zero Problem Time 

Zero problem time is 0500 on the first day of the 
operation. Thus, HHOUR = 300. All insertion times are in 
minutes after 0500, i.e., INSERT(2) = 140 indicates a SIAF 
patrol will be inserted into subarea 2 at 0720 on the first 
day . 

2 . Problem Area Map 

Figure 4 is a map of the problem area translated into 
the coordinate system described in Section III.B. Subarea 
boundaries used are not necessarily exactly as they would be 
in the actual operation, i.e., the river would probably be 
the boundary of patrol subareas on either side of it. The 
following input information is shown on the map for each sub- 
area: 



NR 


XCFF(^IR) 


INSERT (NR) 


XPROBl(NR) 

XPR0B2(NR) 

XPR0B3(NR) 

XPR0B4(NR) 



The values used for expected number of calls for 
fire and target type probabilities are based on the expected 
locations, density, and type of enemy targets as indicated 
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Figure 4. Problem Area Map for STDF Feasibility Study. 




1 




m 



in Section V.B.2. XCFF and XPROBl + XPROB2 are high near the 
river, and XPR0B2 > XPROBl since it is believed that the in- 
surgent force operates primarily in small guerilla groups. 
XPR0B3 + XPR0B4 is high in the mountainous areas to the 
southwest and northwest, indicating suspected movement of 
supplies through those areas, with XPR0B3 > XPR0B4 indicat- 
ing primarily heavy equipment being moved. 

3 . Expected Patrol Durations 

The maximum and minimum patrol durations were 



chosen such that the expected time of extraction of each 
patrol should be approximately 7 days after zero problem 
time. However, patrol duration parameters are given a wider 
range of values in the areas of heavy density of enemy per- 
sonnel. This reflects a higher probability of early extrac- 
tion due to SIAF personnel casualties. 

4 . Assignment Doctrine and Queue Discipline 

The discussion in Section V.B.3.C implies that "any 
ship with guns of adequate range may fire” any mission (ASNDOC 
= 1) , and "personnel missions have priority” in the queue 



(QDISC = 2) . 

5. Refueling Schedule 

The refueling schedule as discussed in Section V.B. 
3.d translates into the following values of refueling sched- 
ule variables: 



REPPOL = 2 

REFSKD = 1 (Computer calculated schedule used) 

COMBO = 60 (Time for round trip between AOE and 
waiting pattern is one hour) 
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NOMO = 720 CNo more refuelings within 12 hours of 
latest extraction time) 

XCRFST(l) =5.25 (5 days, 6 hours after zero prob- 

lem time is 1100 on 6th day) 

XCRFST(2) thru XCRFST(5) = 99.00 

XDDFST(l) =4.25 

XDDFSTC2) = 8.25 

XDDFSTC3) thru XDDFST(5) = 99.00 
CRINT = 240 
DDINT = 240 

Other values of refueling and rearming parameters 
are translated as indicated in Appendix C. 

6 . Firing Time Parameters 

The components of mean firing time as indicated in 
Appendix C were determined as follows: 

a. Denominator of Range Factor 

The range to the target generally has less ef- 
fect on the accuracy of 8''/55 and 6”/47 shells than 5" shells 
since they are larger and designed for longer ranges. 

b. Mean Firing Time Due to Target Type 

A 5"/54 gun is considerably more effective 

against personnel targets and should require considerably 

» 

less time to fire a personnel target mission due to its 
rapid rate of fire. 8"/55 and 6"/47 guns require less time 
against material targets due to their high explosive power, 
while a 5"/54 gun requires less time than a 5”/38 due to its 
rapid rate of fire. 
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D. INVESTIGATION OF VARIOUS NUMBERS AND MIXES OF SHIPS 
PRESENT 

1 . Two Ships 

Rough calculations based on expected number of mis- 
sions, expected service time, and expected problem duration 
indicate that it is very unlikely that any combination of 
two ships would be able to provide adequate support. These 
rough calculations are as follows: 

34 

XCFF(NR) = 414.0 

NR=1 

Expected total service time per mission = 45 minutes 

Expected problem duration = 7 days 

414 missions x 3/4 hours/mission = 310 ship hours 
required 

Two ships X 7 days x 24 hours/ day = 336 ship hours 
available . 

The variability of mission interarrival times, mis- 
sion firing times , and travel distances , and the requirements 
for refueling and rearming would most likely cause an unac- 
ceptable queueing situation to develop. Therefore, combina- 
tions of two ships were not investigated initially. 

2 . Three Ships 

It is not at all apparent from rough calculations 
whether or not three ships would be able to adequately sup- 
port the hypothesized operation. Herein lies the basic 
justification for the use of simulation, since simulation 
accounts for the variability in problem parameters and pro- 
duces an average (expected) result on which to base decisions. 




86 



To investigate the queueing results for three ships 
present, the problem was run 16 times, representing all pos- 
sible combinations of 0 to three ships of each type which 
add to 3, with the requirement that at least one 8'755 ship 
or one 5”/54 ship be present, since the northern extremity 
of the objective area is 25000 yards from the gunline. Ten 
iterations per run (NRIT = 10) were used to determine aver- 
aged results and standard deviations. Table I, on the next 
page, presents the two decision criteria output statistics 
for each of the 16 combinations of ships. Figures in the 
table are the upper and lower limits of a 95% confidence 
interval. These limits were calculated using the equation 

X - 2.26 a= < y < x + 2.26 a= 

X “ ” X 

which was derived as discussed in Section III. I. 4. 

Inspection of the figures in Table I shows that re- 
gardless of the mix of ships, overall ship utilization is 
approximately 50% when three ships are present. From this 
observation it can be readily deduced that utilization when 
four ships are present would be considerably less (probably 
about 35%) while utilization when only two ships are present 
would be greater (probably about 65%) . Based on this con- 
clusion, the simulation was exercised for the seven possible 
combinations of two ships to see if any combination would 
yield a satisfactory average waiting time. The results of 
these seven runs are presented in Table II, on the next 
page . 
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0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

2 

2 

2 

3 

B 

B 

0 

0 

0 

1 

1 

1 

2 



. Decision 
of Three 


Criteria 

Ships 


Output Statistics 


for Combinations 


NR64 7 


NR5 54 


NR5 38 


AVG WAITING TIME 
MIN MAX 


OVERALL SHIP 
UTILIZATION 
MIN MAX 


0 


1 


2 


24.0 


35.0 


.482 


.528 


0 


2 


1 


7.0 


10.2 


.484 


.510 


0 


3 


0 


3.7 


5.3 


.447 


.485 


1 


1 


1 


8.6 


13.4 


.470 


.522 


1 


2 


0 


4.0 


6.2 


.434 


.494 


2 


1 


0 


5.6 


7.8 


.477 


.493 


0 


0 


2 


30.5 


38.1 


.487 


.539 


0 


1 


1 


7.6 


10.2 


.476 


.516 


0 


2 


0 


3.6 


6.6 


.460 


.508 


1 


0 


1 


7.8 


11.2 


.478 


.524 


1 


1 


0 


3.7 


5.7 


.459 


.481 


2 


0 


0 


6.4 


8.0 


.475 


.507 


0 


0 


1 


8.5 


11.1 


.494 


.526 


0 


1 


0 


4.6 


8.4 


.478 


.518 


1 


0 


0 


4.7 


7.3 


.476 


.522 


0 


0 


0 


4.5 


9.9 


.468 


. 508 



I . 


Decision Criteria 
of Two Ships 


Output 


Statistics 


for Combinations 


NR64 7 


NR554 


NR5 38 


AVG WAITING TIME 
MIN MAX 


OVERALL SHIP 
UTILIZATION 
MIN MAX 


0 




1 


1 


48.8 

% 


. 73.4 


.687 


. 743 


0 




2 


0 


26.0 


41.6 


.676 


.742 


1 




1 


0 


29.8 


42.8 


.679 


.715 


0 




0 


1 


65.5 


85.9 


.719 


.767 


0 




1 


0 


25.6 


44.4 


.666 


.732 


1 




0 


0 


33.4 


45.8 


.674 


. 743 


0 




0 


0 


33.1 


46.9 


.685 


.751 
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At this point in the analysis, it was observed that 
no combination of ships present would simultaneously satisfy 
the two decision criteria within the limits of a 95% con- 
fidence interval. IVhen three ships are present, there is no 
combination for which it can be said with 95% confidence that 
OVERALL SHIP UTILIZATION > .50, while there is no combination 
of two ships for which it can be said with 95% confidence 
that AVERAGE WAITING TIME < 10.0 minutes. Thus, one of the 
two decision criteria had to be relaxed. Inspection of 
Tables I and II made it readily obvious that the ship uti- 
lization criteria was the logical criteria to relax. Thus, 
the decision criteria was changed to OVERALL SHIP UTILIZA- 
TION > .43. 



With these revised decision criteria 
nine combinations of three ships satisfy both 
taneously : 



the following 
criteria simul- 



NR855 


NR64 7 


NR554 


NR538 


0 


0 


3 


0 


0 


1 


2 


0 


0 


2 


1 


0 


1 


0 


2 


0 


1 


1 


1 


0 


1 


2 


0 


0 


2 


0 


1 


0 


2 


1 


0 


0 


3 


0 


0 


0 



The striking characteristic of these nine combina- 
tions is that NR538 = 0 in all of them and they are the only 
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nine combinations in which NR538 = 0. This is obviously be- 
cause of the short range of the 5”/38 gun, which limits its 
usability when the objective area is more than 7 miles deep. 

In an attempt to further distinguish between these 
nine combinations of three ships, t\>ro other output statistics 
were investigated. These were MAXIMUM WAITING TIME and PER- 
CENT OF MISSIONS REQUIRED TO WAIT. Table III on the next 
page presents these statistics for the nine combinations of 
ships under consideration. 

Inspection of Table III indicates that two combina- 
tions are considerably less desirable than the others. 

These are; 



NR855 


NR647 


NR554 


0 


2 


1 


1 


2 


0 



The striking characteristic of these combinations 
is that NR647 = 2 in each of them, and they are the only two 
combinations of the original 16 in which NR647 = 2. 

E. CONCLUSION OF THE STUDY 

Based on the foregoing analysis of the output statistics, 
the following conclusion was reached: 

" For a situation such as described herein , any combina- 
tion of three NGFS ships of which none are 5"/38 ships and 
less than two are 6"/47 ships is optimal in terms of providing 
adequate support to SIAF's ashore in an insurgent -contested 
area. " 



90 



Table III. MAXIMUI4 WAITING TIME and PERCENT REQUIRED TO WAIT 
Statistics for the Nine Best Combinations of 
Three Ships 

PERCENT 

MAX WAITING TIME REQUIRED TO WAIT 



NR855 


NR64 7 


NR554 


MIN 


MAX 


MIN 


MAX 


0 


0 


3 


69 


141 


15.5 


32.2 


0 


1 


2 


75 


169 


13.6 


36.6 


0 


2 


1 


183 


251 


21.8 


31.8 


1 


0 


2 


55 


123 


13.7 


39.7 


1 


1 


1 


71 


111 


13.8 


34 . S 


1 


2 


0 


184 


342 


21.6 


33.8 


2 


0 


1 


80 


220 


18.3 


37.7 


2 


1 


0 


88 


174 


18.2 


37.2 


3 


0 


0 


110 


254 


11.7 


41.7 
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In terms o£ STDF feasibility, this can be interpreted 
to mean: 

"If three NGFS ships, of which none are 5"/38 ships and 
less than two are 6"/47 ships, are available for inclusion 
in the STDF, then the STDF concept is feasible in terms of 
providing NGFS to SIAF's ashore in a situation such as 
described herein . " 

It must be emphasized that these conclusions only apply 
to the hypothesized situation. The sensitivity of the con- 
clusions to changes in input parameters has not been tested 
as would certainly be necessary in a detailed feasibility 
study. It is not considered within the scope of this thesis 
to undertake such an extensive study. Rather, it is con- 
sidered that the objectives stated in Section I.C have been 
met at this point. 
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VI. CONCLUSIONS 



It is recognized that there are several ways in which 
the model presented in this paper is capable of refinement 
and extension in order to more accurately simulate realistic 
SIAF support situations and be more useful in investigating 
the overall feasibility of the STDF concept. Although such 
refinement and extensions were not considered to be within 
the scope of this thesis, some of them are pointed out in 
this section for possible future consideration by users of 
the model . 

A. REFINEMENT OF THE MODEL 

1 . Travel Distance Cal culat ions 

The calculations of travel distances on which ship 
assignments are based can be made more realistic by develop- 
ing a dynamic scheme of keeping track of actual positions of 
unassigned ships throughout the run of the model. 

2 . Ship Assignment Doctrines 

Ship assignment doctrines other than the two used 
in this model can be devised to simulate a wider range of 
realistic options available to the STDF commander in as- 
signing ships to missions. 

3. Queue Discipline 

A preempt option for "urgent'' missions can be in- 
troduced to add realism to the model. 

4 . Refueling and Rearming Off-Station Times 

The amount of time spent by each ship off-station 
for rearming and/or refueling can be made more realistic by 



developing a method o£ simulating movement of oilers and am- 
munition ships throughout the problem. 

5 . SIAF Patrol Movements 

The model can be made more realistic by developing 
a dynamic scheme of simulating movement of SIAF patrols with- 
in their assigned subareas. Such a scheme could possibly 
be the linking of this model with the TRW SIAF model in some 
way. 

6 . Distributions of Variable Quantities 

All of the distributions of variable quantities 
used in the model are based primarily on personal experience 
in Naval operations. It is possible that some of the dis- 
tributions do not reflect realistic situations as accurately 
as some other distribution might. Refinement of the model 
should include an investigation of historical or experimental 
data to determine the actual distributions of the variable 
quantities. For example, data concerning actual mission 
firing times should be analyzed to determine if firing times 
actually conform to the exponential distribution or to some 
other probability distribution. 

B. EXTENSIONS OF THE MODEL 

1 . Inclusion of Air Strikes 

The model can be extended to simulate air strikes by 
helicopters or fixed-wing aircraft originating from carriers 
assigned to the STDF when targets could be more effectively 
encountered by air strikes or are beyond the maximum range of 
Naval gunfire. 
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2 . Inclusion of Battleships 

Although there are no battleships in commission at 
the present time, the possibility that battleships could be 
brought back into use in the future, as was the case of the 
NEW JERSEY in the Vietnam War, indicates the desirability of 
including battleships in the model. However, if no ships of 
a particular type were present in a particular problem, the 
model in its present state could be used to simulate the 
presence of battleships by assigning battleship parameters 
to the input variables normally used for ships of the type 
not present. This might require some modification to the 
ship assignment doctrine. 

3 . Multiple Patrols in a Subarea 

The model could be extended to include the option 
of inserting a second patrol into a subarea when the first 
patrol is extracted, and so on. 

C. OVERALL CONCLUSION 

While it is recognized that the model presented in this 
paper is capable of refinement and extension in several ways, 
it is considered that the model in its present state is use- 
ful as a means of providing general ideas concerning the 
queueing aspects of a wide variety of user-defined NGFS/SIAF 
situations. It is also considered that the model can be 
useful in operational planning, i.e., the model could be used 
to investigate the effects of locating the waiting pattern in 
various areas, etc. 
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APPENDIX A: COMPUTER PROGRAM GENERALIZED FLOW CHART 
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I CALCULATE TARGET COORDINATES AND TARGET TYPE 

N ^ y ^ 




DETERXINS MISSION 
TYPE BASED ON TARGET 
RANGE ONLY 



DETERillNE MISSION 
TYPE BASED ON RANGE 
AND TARGET TYPE 



CHECK AVAILABILITY OF SHIPS ELIGIBLE 
TO FIRE THIS TYPE OF MISSION 




CALCULATE TRAVEL 
DISTANCE OF EACH 
AVAILABLE SHIP, 
ASSIGN CLOSEST SHIP 
TO FIRE THE MISSION, 
SET ASSIGNED SHIP 
IN A FIRING STATUS 



I 



CALCULATE SERVICE 
TIME A:^D MISSION 
COMPLETION TIME, 
AND PLACE MISSION 
COMPLETION ON THE 
EVENT CHAIN 



PUT MISSION 
IN THE QUEUE 
AHEAD OF ALL 
OTHERS WITH 
TARGET TYPES 
3 OR 4 



PUT MISSION 
AT END 0 
THE QUEU 
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ADD Dr:rJ2LI2;a TIiV.2 
TO TIME OP COMPIETIOD 

CP rea:-o.:ing, plate 

TOI^TLETICK 07 
TOMB ID ED REFUEL/ 

REARM c:t the e7e:;t 
OHAIN, AIID R-L:0VE 
H.LvRMING CCI.I^'LETIOK 
FROM THE EVERT CHAIN. 
GET CHIP IN A 
" RE7UELIN 3/ REAR!/.ING" 

C T ATITS^ 



C iiXi o U ij A j. -I- 
COMPLETION OF 
REFUELING AND 
PLACE ON THE 
EVENT CHAIN. GET 
SHIP IN A 
"REFUELING" 
STATUS. 



SET SHIP IN A 
"REFUEL AFTER 
MISSION" STATUS 



■© 
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APPENDIX B: COMPLETE COMPUTER PROGRAM LISTING 
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T NUMBER OF MISSIONS AND NUMBER OF EVENTS GENERATED THUS FAR 
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SUBROUTINE SERVIC( SHI P , TAR TYP ,RMG , 1 1 »DIS) 
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SUBROUTINE AVA IL (NRSHI P,SHPTYP ) 
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APPENDIX C: 



PRINTOUT OF INPUT DATA FOR STDF FEASIBILITY STUDY 
(For the Case When NR855 = NR554 = 1, NR558 = 0) 



INPUT 



NUMBER 
NUMBFP 0*= 
NUMBFR OF 
NUMBER 0^ 



8«/55 SHIPS PRESENT 
6**/47 SHIPS PRESENT 
5«/54 SHIPS PRESENT 
5«/38 SHIPS PRESENT 



1 

1 

1 

0 



MAX RANGE OF 8»/55 GUN = 30000 YDS 
MAX RANGE OF 6”/47 GUM = 21000 YDS 
MAX RANGE OF 5*V54 GUN = 26000 YDS 
MAX RANGE CF 5”/38 GUN = 15000 YDS 



DISTANCE OF WESTERN EDGE OF WAITING PATTERN FROM Y-AXIS = 10000 YDS 
DISTANCE OF NORTHERN EDGE OF WAITING PATTERN FROM GUNLINE = 5000 YDS 

E-W LENGTH OF WAITING PATTERN = 20000 YDS 
N-S LENGTH OF WAITING PATTERN = 5000 YDS 



SPEED TO FIRING POSITION = 20 KNOTS 
MIN SETUP TIME = 3 MINS 

MAX SETUP time = 15 MINS 



ASSIGNMENT DOCTRINE: 1 

QUEUE DISCIPLINE: 2 

REPLENISHMENT POLICY: 2 

TYPE OF FUELING SCHEDULE USED: 2 



♦ COMPONENTS OF MEAN FIRING TIME « 

MEAN FIRING TIME DUE TO TARGET TYPE: (MINUTES FPCM TABLE BELOW) 



TARGET TYPE 



SHIP 


1 


2 


3 


4 


1 


15 


10 


12 


10 


2 


14 


8 


11 


9 


3 


8 


5 


14 


10 


4 


18 


12 


18 


12 



DENOMINATOR GF RANGE FACTOR: (FROM TABLE BELOW) 



TYPE SHIP 


1 


TARGET 

2 


TYPE 

3 


4 


1 


10000 


8000 


6000 


6000 


2 


lOOCC 


8000 


6000 


6000 


3 


10000 


8000 


6000 


6000 


4 


8000 


4000 


4000 


4000 
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AREA 

NR 

I 

2 

3 

A 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 



AREA INFGRMaTION 



EXP PATROL TARGET TYPE 





AREA eOUMDARIBS 


DURATION 


EXP NR 


PROBABILITIES 


NORTH 


SOUTH 


WEST 


EAST 


MIN 


MAX 


OF CPF 


1 


2 


3 


4 


25000 


20000 


0 


4000 


6.0 


7.0 


8.0 


0.10 


0.20 


0. 50 


0.20 


250GC 


2CC00 


4000 


8000 


6 .4 


7.4 


8 .0 


0.10 


0.20 


0.50 


0.20 


25000 


2CCOO 


8000 


12000 


6.0 


7.0 


8. 0 


0. 10 


0.20 


0.50 


0.20 


250CC 


20000 


12000 


16C00 


5.8 


6.8 


8.0 


0.20 


0.30 


C.4C 


0. 10 


25000 


2C000 


16000 


19000 


4.6 


7.6 


10.0 


0.20 


0.40 


0.10 


0.30 


25000 


20000 


19000 


22000 


5.0 


7.0 


1C. 0 


0.20 


C.4C 


0. 10 


C.30 


2500C 


2C000 


22000 


28000 


5.9 


6.9 


12.0 


0.20 


0.20 


0.50 


0.10 


25000 


20000 


28000 


40000 


6.5 


7.0 


12.0 


0.10 


0.15 


0.65 


0.10 


20000 


15000 


0 


4000 


6.5 


7.5 


11.0 


O.IO 


0.20 


0.50 


0.20 


200C0 


15000 


4000 


8000 


5 .7 


6.7 


12.0 


0.10 


C .20 


0.50 


0.20 


2000C 


15000 


8000 


13000 


6.2 


7.2 


15.0 


0.10 


0.10 


0.60 


0.20 


20000 


15000 


13000 


18000 


5.0 


6.5 


18.0 


0.25 


0.25 


C.35 


C. 15 


20000 


15000 


18000 


21000 


4.5 


7.5 


15.0 


0.20 


0.40 


0.10 


0 .30 


20000 


15000 


21000 


2400C 


5 .3 


7.3 


15.0 


C.20 


0.40 


C. 10 


0.30 


20000 


15000 


24000 


3200C 


6.5 


7.5 


15.0 


0.10 


0.15 


0.40 


0.35 


2 0000 


12000 


32000 


40C00 


6.6 


6. 8 


10.0 


0.10 


0.20 


0.10 


0.60 


15000 


10000 


0 


4000 


6.5 


7.5 


11.0 


0. 10 


0.20 


C. 50 


0. 20 


1500C 


ICOOC 


4000 


8000 


6 .0 


7.0 


12.0 


0.15 


0.25 


0.30 


0.30 


15000 


12000 


8000 


13000 


5.8 


6. 8 


10.0 


0. 10 


0. 10 


0.60 


0.20 


150C0 


12000 


13000 


18000 


5.7 


6.7 


11.0 


0.10 


0.3C 


0.50 


0. 10 


150C0 


12000 


18000 


22000 


4.5 


6.5 


18.0 


0.20 


0.40 


0.10 


0.30 


15000 


12000 


22000 


26000 


5.0 


6. 8 


15.0 


0.20 


C.4C 


0. 10 


0.30 


15000 


12000 


260C0 


32000 


6.0 


7.0 


12.0 


0.20 


0.3C 


0.2 5 


0.25 


10000 


6000 


0 


8000 


6.0 


7.0 


8.0 


0.05 


0.35 


0.20 


0.40 


12000 


8000 


8000 


15000 


6.2 


7.2 


IC.O 


0.10 


0.25 


0.35 


0. 30 


12000 


8000 


15000 


21000 


5.5 


7.0 


11.0 


0.20 


0.40 


0.20 


0.20 


12000 


8000 


210C0 


24000 


4.0 


7.0 


22. 0 


0.20 


C.4C 


0. 10 


0.30 


120C0 


8000 


24000 


28000 


4.5 


7.5 


18.0 


0.20 


0.40 


C. 10 


0.30 


12000 


8000 


28000 


4C000 


6.4 


6.6 


7.0 


0.10 


0.20 


0.10 


0.60 


8000 


4000 


8000 


15000 


6.2 


7.2 


7. 0 


0.05 


0. 35 


0.20 


0.40 


8000 


4000 


15000 


21000 


5.0 


6.5 


8.0 


0.20 


0.50 


0.10 


0.20 


8000 


5000 


21000 


26C00 


4.0 


7.5 


18.0 


0.40 


0.30 


0.20 


0. 10 


8000 


5000 


26000 


30C00 


4.0 


7.0 


23.0 


0.40 


0.30 


0. 20 


0. 10 


80C0 


5000 


30000 


40000 


6,3 


6.5 


6 .0 


0.10 


0 .20 


0.10 


0.60 
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♦ REFUELING AND REARMING DATA * 



GPF-STATION TIME PARAMETERS: 

(MINUTES Of=F STATION PER EVOLUTION FROM TABLE BELOW) 



CRUI SER S 

REFUEL REARM 

MIN MAX MIN MAX 


DESTROYERS 
REFUEL REARM 

MIN MAX MIN MAX 


150 180 


240 300 


120 150 


180 240 


INITIAL ammunition LEVELS 
(ROUNDS ON BOARD FROM 


’ BELOW T ABLE) 




SHIP NR 


8''/55 6» 


SHIP TYPE 
/47 5"/54 


5"/ 38 


1 


8000 10000 12000 


6000 



MINI MUM 


ALLOWED AMMUNITION LEVELS = 


8»/55 : 


1000 RDS 


6"/47 : 


1000 RDS 


5”/54: 


2000 RDS 


5»/38: 


1000 RDS 



EXPECTED AMMUNITION 


EXPEND 


ITURFS: 








(ROUNDS 


PER 


MINUTE PER 


MI SSION) 
















TARGET 


TYPE 






TYPE SHIP 




1 




2 


3 








MIN 


MAX 


MIN 


MAX 


MIN M 


AX 


MIN 


1 


4 


10 


3 


9 


2 


8 


2 


2 


4 


11 


3 


10 


2 


8 


2 


3 


6 


22 


5 


20 


5 


18 


5 


4 


4 


10 


3 


9 


3 


10 


3 
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APPENDIX D: PRINTOUT OF MISSION HISTORY FOR STDF FEASIBILITY STUDY - (FIRST 

ITERATION WITH INPUT PARAMETERS lAW APPENDIX C) 
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AVERAGE TRAVEL TIME = 16.3 MINS 

AVERAGE FIRING TIME = 11.9 MINS 

average total SERVICE TIME (TRAVEL + SETUP FIRING TIME) 
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OVERALL SHIP UTILIZATION 
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APPENDIX F; PRINTOUT OF OVERALL STATISTICAL RESULTS FOR STDF 
FEASIBILITY STUDY (WITH INPUT PARAMETERS lAW 
APPENDIX C) 



****** OVERALL STATISTICAL RESULTS 






FOR 10 



_ _ OF 

iterations 



THE MEAN 



QUANTITY 

TOTAL PROBLEM TIME (DAYS) 

TOTAL NUMBER CF MISSIONS FIRED 



MEAN 


STD DEV 


7.34 


0.04 


385 


7 



AVERAGE WAITING TIME PEP MISSION 
missions which WERE REQUIRED 

AVERAGE WAITING TIME PER MISSION 

MAXIMUM WAITING TIME (MINS) 



(INCL ONLY 
TO WAIT) (MINS) 

(ALL MISSIONS) 



19.9 



5.0 

95 



1.78 

0. 71 
14 



AVERAGE NUMBER OF MISSIONS AWAITING ASSIGNMENT 
MAXIMUM NUMBER OF MISSIONS AWAITING ASSIGNMENT 

peouieed 

’’^"'^"FOR^SsiGNHilSr FEOUIREO 10 UAIT 



0.18 


C.03 


5 


0.54 


95 


6.71 


24.4 


4. 55 



ONLY 

'‘’'"'‘(IePSLhonsY' 

MAXIMUM TIME FROM CFF UNTIL FIRST SHOT 



45.4 

25.7 

30.6 

132 



2.12 
0.07 
0.78 
13. 5 



AVERAGE TRAVEL TIME (MINS) 

AVERAGE FIR ING TIME 

AVERAGE TOTAL SERVICE TIME 

(TRAVEL + SETUP + FIRING TIME) 



16.6 C. 13 
11.8 0.27 
37.4 0.35 
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OVERALL SHIP LTILIZATION : MEAN = O.A69, STO DEV - 0.007 
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